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ABSTRACT 


Stellar X-ray and UV radiation can significantly affect the survival, composition, and long-term 
evolution of the atmospheres of planets in or near their host star's habitable zone (HZ). Especially 
interesting are planetary systems in the solar neighborhood that may host temperate and potentially 
habitable surface conditions, which may be analyzed by future ground and space-based direct-imaging 
surveys for signatures of habitability and life. To advance our understanding of the radiation environ- 
ment in these systems, we leverage ~3 Msec of XMM-Newton and Chandra observations in order to 
measure three fundamental stellar properties at X-ray energies for 57 nearby FGKM stellar systems: 
the shape of the stellar X-ray spectrum, the luminosity, and the timescales over which the stars vary 
(e.g., due to flares). These systems possess HZs that will be directly imageable to next-generation 
telescopes such as the Habitable Worlds Observatory and ground-based Extremely Large Telescopes 
(ELTs). We identify 29 stellar systems with Lx/Lyo| ratios similar to (or less than) that of the Sun; 
any potential planets in the habitable zones of these stars therefore reside in present day X-ray ra- 
diation environments similar to (or less hostile than) modern Earth, though a broader set of these 
targets could host habitable planets. An additional 19 stellar systems have been observed with the 
Swift X-ray Telescope; in total, only ~30% of potential direct imaging target stars has been observed 
with XMM-Newton, Chandra, or Swift. The data products from this work (X-ray light curves and 
spectra) are available via a public Zenodo repository (doi: 10.5281/zenodo.11490574). 
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1. INTRODUCTION 


In the coming decades, ground- and space-based di- 
rect imaging will afford the best potential for detecting 
and spectrally characterizing terrestrial planets in the 
habitable zones of our nearest neighbors (Turnbull et al. 
2021; Kopparapu et al. 2018; Li et al. 2021; Kane et al. 
2018; Currie et al. 2023). While the James Webb Space 
Telescope (JWST) has and will continue to characterize 
some rocky exoplanets (e.g., Greene et al. 2023; Zieba 
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et al. 2023; Lincowski et al. 2023; Lustig-Yaeger et al. 
2023; Moran et al. 2023), its ability to spectrally ex- 
amine temperate terrestrial—and therefore potentially 
habitable—exoplanets is limited primarily to target sys- 
tems with transiting planets and late M dwarf host stars. 
This is due to the lack of coronagraphic instrumentation 
on JWST to image rocky planets around more massive, 
Sun-like stars, and because the characterization of ter- 
restrial planets through transit observations is most fa- 
vorable for M dwarf hosts, which have a combination 
of favorable transit probabilities, planet-star size ratios, 
and planet-star separations that allow for the probing 
of low atmospheric altitudes when accounting for the ef- 
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fects of refraction (García Muñoz et al. 2012; Bétrémieux 
& Kaltenegger 2013, 2014). Moreover, some key candi- 
date biosignature molecules like O2 cannot be detected 
directly at Earth-like abundances with JWST via transit 
spectroscopy even for the most favorable targets (e.g., 
Lustig-Yaeger et al. 2019; Pidhorodetska et al. 2020; 
Meadows et al. 2023), but could be more easily detected 
by reflected light observations of planets in nearby sys- 
tems (The LUVOIR Team 2019; Gaudi et al. 2020). 

The occurrence rates derived from Kepler indicate 
that terrestrial planets are likely present within the hab- 
itable zones (HZs) of a large fraction of nearby stars 
(e.g., Bryson et al. 2021; Kunimoto & Matthews 2020; 
Dattilo et al. 2023; Fulton et al. 2017; Burke et al. 2015; 
Kane et al. 2014). The Astro2020 Decadal Survey has 
identified the discovery and characterization of these 
worlds as the highest priority for NASA’s investments 
(National Academies of Sciences & Medicine 2023), and 
the Habitable Worlds Observatory will build on the LU- 
VOIR and Habitable Exoplanets Observatory concepts 
(The LUVOIR Team 2019; Gaudi et al. 2020) to real- 
ize this goal. Due to detection biases, the vast majority 
of nearby temperate rocky planets amenable for direct 
imaging have not yet been discovered, but will be re- 
vealed by HWO in its survey phase (e.g., The LUVOIR 
Team 2019; Gaudi et al. 2020), precursor science with 
extreme precision radial velocity (EPRV) measurements 
(Crass et al. 2021; Morgan et al. 2021), or direct imaging 
with upcoming 30-m class telescopes (Fujii et al. 2018; 
Currie et al. 2023). 

A detailed understanding of the target stars them- 
selves is a strong prerequisite for maximizing the science 
yield of exoplanet imaging missions, including extracting 
the statistical features of the overall planet population 
(e.g., Van Eylen et al. 2014) as well as characteriz- 
ing the climates, atmospheric composition and stability 
(Linsky & Redfield 2024; Louca et al. 2023), and po- 
tential biosignatures of individual planets (Rugheimer 
et al. 2015; Schwieterman et al. 2018). Evaluating plan- 
etary habitability, specifically, must consider myriad fac- 
tors beyond circumstellar location including the planet's 
bulk composition (Hinkel & Unterborn 2018; Maroun- 
ina & Rogers 2020; Unterborn et al. 2023) and stellar 
X-ray and extreme ultraviolet (XUV; 1-912 À; corre- 
sponding to an energy range of 0.01-12 keV) emission 
(France et al. 2016; Garcia-Sage et al. 2017; Luger & 
Barnes 2015; Richey-Yowell et al. 2019). Planets or- 
biting in or near their star's HZ can experience high 
levels of XUV radiation, which directly influences atmo- 
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spheric survival, composition, and long-term evolution 
(Peacock et al. 2020; Johnstone et al. 2021). In some 
cases, planetary atmospheres and surface volatiles can 
be entirely driven away by XUV irradiation (e.g., Sanz- 
Forcada et al. 2011; Fromont et al. 2024; Luger & Barnes 
2015). Further, the atmospheric states of planets orbit- 
ing M dwarf stars may be profoundly different than for 
those orbiting Sun-like FGK stars due their to extended 
super-luminous pre-main sequence phase (Baraffe et al. 
2015) during which XUV radiation could ablate away 
planetary atmospheres or lead to atmospheric oxidation 
via massive H escape (Garcia-Sage et al. 2017; Luger & 
Barnes 2015; Krissansen-Totton & Fortney 2022). Out- 
gassing could replenish these atmospheres (Kite & Bar- 
nett 2020; Swain et al. 2021; Unterborn et al. 2022; Kane 
et al. 2020a) though the retention and photochemistry 
of the atmosphere would continue to be controlled by 
the host star spectrum and activity (Roettenbacher & 
Kane 2017; France et al. 2020, 2016; Loyd et al. 2018a,b; 
Do Amaral et al. 2022). 

While X-ray detections have been made for some 
HWO target stars (Harada et al. 2024) and more de- 
tailed studies of the X-ray properties of low-mass nearby 
stars have been performed (Brown et al. 2023; France 
et al. 2020; Loyd et al. 2018b), no comprehensive and 
uniform X-ray analysis of potential direct imaging target 
stars has yet been done. In this paper, we present a com- 
prehensive analysis of archival X-ray observations for all 
FGKM stellar systems with HZs that will be directly- 
imageable with upcoming telescopes such as HWO and 
ground-based Extremely Large Telescopes (ELTs). A 
follow-up paper yielding panchromatic XUV-IR spectra 
(including modeled EUV wavelengths) for each one of 
these systems is forthcoming (Peacock et al., in prep). 
We describe our sample and the archival data in Sec- 
tions 2 and 3. Our analysis procedure is detailed in 
Section 4. We extract light curves for stars detected at 
high significance and identify periods of count rate vari- 
ability (such as flaring events or epochs of elevated count 
rates). We then extract X-ray spectra for stars, sepa- 
rating variable- from non-variable epochs, and fit the 
resulting spectra with one-, two-, or three-component 
thermal plasma models. Finally, the average best-fit 
spectral models are used to convert observed count rates 
(or count rate upper limits, in the case of non-X-ray- 
detected stars) to luminosities. In Section 5, we discuss 
the evolution of the Lx / Lo) ratio with effective temper- 
ature and stellar age, and present an extensive appendix 
describing the X-ray observations and a summary of the 
known stellar physical parameters and exoplanet sys- 
tems from the literature for each stellar system in our 
sample. 


2. THE SAMPLE 


To identify the highest priority targets for “deep 
dive” investigations of nearby stars and their plane- 
tary systems, we have used the predicted imaging and 
spectroscopy capability of (1) NASA’s next generation 
space-based Great Observatories such as HWO, and (2) 
ground-based ELTs. To first order, the performance of 
these facilities is described in terms of two specifications: 
the inner working angle (IWA) — the smallest planet- 
star angular separation that can be imaged (expressed 
in milliarcseconds, or mas), and the minimum planet to 
star flux ratio (or fractional planet brightness, FPB) — 
where starlight suppression technology and postprocess- 
ing algorithms can separate the planet signal from stellar 
photons and systematic noise. Intriguingly (and some- 
what vexingly), for HZ planets the IWA requirement 
favors stars that are more luminous and therefore have 
wider HZs, while the limiting flux ratio favors target 
stars that are less luminous and therefore outshine their 
HZ planets by a smaller factor (Turnbull et al. 2012). A 
third specification, planet limiting magnitude (Vim), de- 
termines which of these targets planets would be bright 
enough to obtain spectra in a reasonable amount of ob- 
serving time — and this requirement simply favors the 
nearest stars. 

To form our preliminary target list, for space-based 
missions we chose stars for which the HZ (0.95 - 1.7 
AU scaled by Ls) falls at least partially outside an 
inner working angle of 58 mas, and for which Earth- 
sized planets of albedo —0.3 would be brighter than 
Vim < 30, with a fractional planet brightness (FPB) 
>4x 10711, For ground-based observatories, we exam- 
ined the 25 stars within 10 pc having the most favor- 
able FPB whose HZs are at least partially visible out- 
side an IWA of ~30 mas, consistent with the literature 
(Brandl et al. 2018; Quanz et al. 2015). This results 
in a well-defined list of 175 target stars whose HZs fall 
at least partially within the detection space of upcom- 
ing exoplanet direct imaging observatories. We com- 
pared our target list with the potential target stars for 
HWO (Mamajek & Stapelfeldt 2024), who used slightly 
different cut-off requirements for Viim, FPB and IWA. 
Our list contains more M-dwarfs (more favorable targets 
for ELTs) while the Mamajek & Stapelfeldt (2024) list 
includes more F stars (more favorable for space-based 
observations). Mamajek & Stapelfeldt (2024) also ex- 
cludes stars in close binaries or higher-order multiple 
systems, while we produce a more inclusive X-ray cata- 
log of nearby stars, assuming technical challenges posed 
by binaries may eventually be surmounted. There are 
229 unique stars when combining the two lists, which 
we summarize in Table 1. 
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Table 1. Potential Target Stars for Future HZ Imaging Surveys 


HD Common Distance Spectral From Potential 


ID Name (pc) Type This Work? HWO Target® 


166 V439 And 13.77 G8V "4 "4 
693 6 Cet 18.89 F8V "4 
739 | 0 Sd 21.72 F5V "4 
1326 GX And 3.57  MI1.5V "4 

F9.5V "4 v 


1581 ¢ Tuc 8.61 


NoTE—^From Mamajek & Stapelfeldt (2024). Select entries 


are shown to illustrate the table form and content. The full 


machine-readable table (229 rows) is available online from the 
journal. 


Using this input list of 229 viable targets for direct 
imaging campaigns, we searched the archives of the 
XMM-Newton Observatory and the Chandra X-ray Ob- 
servatory for publicly available data. Many M stars (and 
some K stars) in this initial target list were previously 
studied as part of the Measurements of the Ultraviolet 
Spectral Characteristics of Low-mass Exoplanetary sys- 
tems (MUSCLES) and Mega-MUSCLES Hubble Space 
Telescope Treasury programs, which characterized the 
X-ray/UV emission from K and M exoplanet host stars. 
We incorporate these results, presented in Brown et al. 
(2023), into our work here (see Section 5). We retrieve 
192 XMM-Newton and Chandra observations for 57 stel- 
lar systems containing Sun-like and low-mass stars from 
the X-ray archives, amounting to ~3 Msec of observing 
time (not including the MUSCLES/Mega-MUSCLES 
observations). An additional 19 stellar systems have 
been observed by Swift. Figure 1 shows a map of the 
3D spatial distribution of these nearby stars and indi- 
cates which stars have X-ray observations available. 

There are three M-dwarfs highlighted in Mamajek & 
Stapelfeldt (2024) as strong potential HWO targets (La- 
caille 8760, Lacaille 9352, and Lalande 21185). Both 
Lacaille 8760 and Lacaille 9352 were observed by Swift 
and are included in this work. A detailed study of the 
panchromatic spectrum and habitability potential of La- 
lande 21185 is in preparation (Schwieterman et al., in 
prep) and is not included in the sample here. Roughly 
half of the target stars contained in our study are F- 
and G-type stars. Table 2 provides a summary of the 
stellar physical parameters for all stars included in our 
study, including the effective temperature (Teg), spec- 
tral type (from Gray et al. 2006), stellar mass (M), 
bolometric luminosity (Lpo1), and approximate age (in 
Gyr). All data were compiled using the NASA Exo- 
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planet Archive”. Right ascension (R.A.), declination 
(Dec.), and distances were taken from Gaia DR3 (Gaia 
Collaboration et al. 2023). More detailed discussion of 
each stellar system is provided in Appendix A. 


Figure 1. Map of nearby potential direct imaging target 
stars, color-coded by spectral type. Stars with blue halos 
have available X-ray observations. The yellow star at the 
center of the diagram represents the Sun, with white con- 
centric rings showing distances of 5 pc, 10 pc, and 15 pc. 
An animated version of this figure (8 seconds in duration) is 
available from the online journal showing the 3D distribution 
through different elevation and azimuthal angles (colors and 
symbols are the same as the static figure). 


3. X-RAY DATA 


We retrieved all publicly available imaging observa- 
tions for our target stars from the XMM-Newton? and 
Chandra* archives. Table 3 summarizes the number of 
observations and the approximate exposure times uti- 
lized in this work. Below, we provide a more detail de- 
scription of how the observations were processed and 
the methods used to extract light curves and spectra for 
each facility. 


2 See https: //exoplanetarchive.ipac.caltech.edu/ 
3 See https://nxsa.esac.esa.int /nxsa-web/#home 
4 See https: //cda.harvard.edu/chaser/ 


5 https: //www.cosmos.esa.int /web/xmm-newton/sas 


6 See the XMM-Newton ABC Guide at https://heasarc.gsfc.nasa. 


gov/docs/xmm/sl/epic/image/sas_cl.html 


3.1. XMM-Newton 


We retrieved XMM-Newton observations from the 
XMM-Newton Science Archive and reprocessed the data 
from the raw evti files using SAS 18°. The Euro- 
pean Photon Imaging Camera (EPIC) on board XMM- 
Newton carries a set of three X-ray CCD arrays: one 
array uses pn-CCDs (referred to as PN) and two uti- 
lize metal oxide semi-conductor CCDs (referred to as 
MOS1 and MOS2). The PN data were reprocessed 
with epproc and the MOS data were reprocessed us- 
ing emproc and filtered on an energy range of 0.2-15 
keV. Due to the brightness of many of the stars, pile- 
up (when two or more photons fall on a pixel in less 
than the readout time) is a significant concern. The 
PATTERN parameter records the number and pattern 
of CCD pixels that are triggered for a given event; we 
use the recommended PATTERN=0 to mitigate pile-up 
effects?. Background light curves were extracted using 
evselect and inspected for background flares, and good 
time intervals (GTIs) were generated to reject observa- 
tion times exhibiting strong background flaring events. 

We inspected the reprocessed, cleaned images for ev- 
idence of an X-ray source at the location of the target 
star. If the target star was detected at 2500 counts in 
the PN image, we extracted source light curves (binned 
to 100 s) to search for count-rate variability. We used 
evselect to extract both PN and MOS1/2 spectra for 
stars with 22000 net counts, rmfgen to generate redis- 
tribution matrix files (RMFs), and arfgen to generate 
ancillary response functions (ARFs) for all available ob- 
servations. All spectra were binned to contain at least 
25 counts per bin. 

The majority (~83%) of the observations analyzed in 
this work are from the XMM-Newton archive, contain- 
ing ~75% of the total exposure time. However, with 
a spatial resolution of ~15”, only stellar binaries with 
angular separations larger than —30" can be resolved 
with the PN camera. Figure 2 (left panel) shows a 10 ks 
XMM-Newton/PN image of the binary system GL 412, 
where the two components (and a third, unidentified X- 
ray source) are resolved. 


Table 2. Properties of Sample Stars 


Star Alternate R.A. Dec. Spectral Distance Mass Tett Lio Age* 
Name Name(s) (32000) Type* (pc) (Mo) (K) (Lo) (Gyr) References 
Q) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 
n Crv HIP 61174 12:32:03.75 -16:11:46.55 F2V 17.96t0.17  1.48*029 68541124 4.62*0-22 15707 1,2,50 
€ Oph HIP 84893 17:21:00.68 -21:06:49.81 F2V 17.0640.10 1.444938 = 67564132 433t017 1,2 
v And A HD 9826 01:36:47.60 41:24:13.54 — F8V 13.4140.06 — 1.290.04 61574112 d3bhols 3.12 1,2,29 
Hor HR 810, HD 17051 02:42:34.03 -50:47:57.54 F9V . 17.32:0.02 1.174917 61464129 1:741 te 0.47; 2.72 1-4,9 
v Phe HIP 5862 01:15:12.13 -45:31:51.14 FOV — 15.17+0.04 1.1570 19 6119 739 2.00+0.06 2.88 1,2,3,9 
y Pav HIP 105858 21:26:26.81 -65:21:45.50 FOV = 9.2740.04 1.151019  6109*119 1.47£0.05 1,7.25 1,3,44,45 
B Vir HIP 57757 11:50:42.51 01:45:48.67 FOV  11.12*005 . 1137018 60711320 grim 2.96 1,2 
GJ 1095 HIP 35136 07:15:50.19 +47:14:20.89 FOV 16.861010 1.12 5866 1.54+0.06 1,7 1,2,10 
LHS 237 GJ 288, HIP 37853 07:45:34.63 -34:09:53.86 F9V 5.20+0.13 1.11 5744 1.51+0.06 1,2 
l Per HIP 14632 03:09:06.10 +49:36:46.35 F9.5V  10.51+0.03 1.0879 15 5938+144 2234008 4,.20+0.48 1,2,9 
B Hyi HIP 2021 00:25:55.80 -77:15:10.06 GOV = 7.4750.03.— 1.14*0 05 5839 3.7040.05 6.32 1,2,9 
B Com HIP 64394 13:11:51.43 +27:52:55.58 GOV  9.18+0.03 — 1.10*01$ 5969+188 115905 «1.2, 1.7 1,2,9 
LHS 208 HD 39091, r Men  05:37:11.89 -80:27:52.05 GOV 8.27+0.02 — 1.094--0.039 6037+45 1.44+0.02 2.08* 120 5 
p CrB GJ 9537, HIP 78459 16:01:02.41--33:18:00.28 GOV 7.4740.02 1.05*012 5833415) 1.817006 11.04 1,2,9 
GL 672 HIP 84862 17:20:39.74 +32:27:47.22 GOV 4.54+0.02 1.054016 58171120 1.297008 12.04 1,2,9 
x! Ori HIP 27913, GJ 222 05:54:22.76+20:16:33.07 GOV  8.66+0.08  1.028*0 05, 5822 1.17+0.06 0.35; 4.32 1,2,6,9 
GL 788 HIP 100017 20:17:32.60 +66:51:18.04 GOV 7.57£0.52  1.0240.03 5886 1.18+0.10 3.767155 2,9 
47 UMa GJ 407, HD 95128  10:59:27.53--40:25:49.80 GOV 3.80--0.03 1.005+0.047 5829+95 1.58+0.04 6.48* 1 04 9,29 
GL 620.1A HD 147513 16:24:01.39 -39:11:34.68 G1V  12.90340.02 1.061317 58734103 1.03278 004 2;<1;0.4 1-4,9,19 
GL 311 x! UMa, HD 72905 08:39:11.64+65:01:16.67 G1.5V — 14.27*0 14 1.01 5920 1.03+0.05 0.2 2,6,8 
a Cen A GJ 559A, LHS 50 — 14:39:36.49 -60:50:02.37 G2V — 1.347x0.003 — 1.06*0 05 5801 1.6140.07 T BATTS 2,9 
44 Boo A Boo, HIP 73695 A 15:03:46.74+47:39:15.89 G2V 12.8+0.2 1.04+0.1 5877 1.552 1.5 46,47,48 
18 Sco GJ 616, HIP 79672 16:15:37.52 -08:22:17.91 G2V 4.1340.02 — 1.04'0 14 5791155" 1.09* 0:05 5.84* 195 1,9,10 
51 Peg HIP 113357, GJ 882 22:57:28.22.-20:46:08.78 G2V 5.4640.03 1.03*0 05 57587193 1:37: 50:02 6.76 1,2 
C! Ret LHS 171, HIP 15330 03:17:49.26 -62:34:20.76 G2V 2.12+0.08 — 0.96*005 5699 0.80-£0.05 1.56£0.44 2,9 
HD 136352 GJ 582, v? Lupi — 15:21:45.55 -48:19:07.88 G2V 4.6823-0.02 0.87+0.04 5664+61 1.04+0.06 123432 11,12 
GL 327 HIP 43726 08:54:17.50 -05:26:03.56 G3V 17.42'0 47  1.002-0.04 5790 1.09*0 05 1.32 2,9 
u Ara HD 16069 17:44:08.68 -51:50:05.65 G3IV/V  15.6040.03 1.034915  5772*11? 1.90* 0:05 5.7£0.6 12/7 
k! Cet HD 20630 03:19:21.98 +03:22:14.22 G5V  — 9.144-0.03 1:02:0:16-. ..5712:5187 0.86* 0-02 0.35; 0.7; 2.2 1,2,6,20,21 
HD 140901 GL 599 15:47:28.54 -37:55:02.15 G7IV/V 15.25+0.01 1.06+0.05 55874153 1.22+0.03 3.242.8 1,49 
GJ 777A LHS 3510, HD 190360 20:03:38.25 +29:53:40.08 G7IV/V — 16.012-0.02 0.990.08 5552+11 1.20+0.06 13.4 1,2,9 
£ Boo A 37 Boo, HD 131156 14:51:23.53+19:06:00.97 G7V — 6.71£0.16  0.9040.04 5551-21 0.55-£0.04 <4 38 
6 Pav GJ 780 20:08:46.81 -66:11:13.53 G8IV  6.09+0.01  0.99*915  5590*128 1.25*0 02 6.2 1,2 
LHS 2156 HIP 47080 09:35:38.55 +35:48:32.33 GSIV/V 11.20+0.02 — 0.96*0!6 5499+156 0.80* 004 1,2 
82 Eri GJ 139, HD 20794  03:20:00.08 -43:03:59.50 G8V  6.00340.009 0.94+0.12 5413191, 0.65+0.01 5.76; 2 12 — 1,2,9,13 
T Ceti HD 10700, GJ 71  01:44:02.17 -15:56:01.25 G8.5V  3.60340.007  0.92*012 53331124 0.51£0.01 >12 1,2,9 
44 Boo Ba 15:03:46.74 +47:39:15.89 KOV 12.8+0.2 0.98 5300 0.51 46,47,48 
55 Cnc A p! Cnc, HIP 43587 08:52:35.22 +28:19:47.22 KOIV/V 12.59+0.01 0.907955 52521373 0.64+0.02 9.5 1,2 
70 Oph A LHS 458, HIP 88601A 18:05:27.47 +02:29:42.81 KOV  5.15x0.03  0.894-0.02 5300+50 0.65+0.02 6.2+1.0 1,22 
40 Eri A GJ 166A, HD 26965 04:15:13.91 -07:40:05.08 K0.5V  5.0440.01 0.85*0 1, 5092+149 0.43£0.0 6.94.7 1,14,23 
6 Eri HIP 17378 03:43:14.80 -09:45:36.30 K1II/IV 9.04+0.07 1.19 5095 g.3r 007 6.28 2 
GL 451A Groombridge 1830 11:53:04.16 +37:41:34.18 K1IV —— 9.16*0.18 0.66 4950 0.23--0.04 5.00.3 2,6 
GL 117 HIP 13402, HD 17925 02:52:32.56 -12:46:14.00 K1.5V 10.355+0.005 0.891993 5225793, 0.40+0.0 1.5;<1.2;0.1  1,2,6,9 
a Cen B GJ 59 B, LHS 51 — 14:39:35.06 -60:50:15.10 K2IV  1.347+0.03 — 0.87-:0.07 5234+63 0.52--0.05 8, >11 2,9,30 
GL 783 A HD 191408A, LHS 48620:11:12.54 -36:06:20.58 K2.5V 6.05+0.03 0.74 4922 0.28+0.10 7.05£0.65 — 1,2,6,25 
GL 183 HD 32147, HIP 23311 05:00:49.59 -05:45:30.96 — K3V  8.845+0.003 0.811955 49317738 0.29+0.0 2 1,3,6 
GL 892 HD 219134 23:13:21.06+57:10:10.80 K3V  6.531£0.004 0.7940.03 4817+62 0.29+0.0 11.0+2.2,12.46 1,2,36 
GJ 667 A HIP 84709 17:18:56.53 -34:50:25.40 K3V 6.97+0.83 0.59 4810 0.13*0-17 2,3,31 
LHS 1875 GL 250, HIP 32984 06:52:17.47 -05:10:25.42  K3.5V — 8.747+0.004 — 0.76*0 09 — 47191137 0.225002 3.1,«0.48 1,2,9 


Table 2 continued 


Table 2 (continued) 
Star Alternate R.A. Dec. Spectral Distance Mass Test Lbol Age* 
Name Name(s) (J2000) Type“ (pc) (Mo) (K) (Lo) (Gyr) References 
Q) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 
GL 570 A HIP 73184 14:57:29.18 -21:25:23.31  K4V — 5.881-0.003 0.754995 46817188 0.22 * 001 «0.6;3 1,2,9 
£ Boo B 14:51:23.19 -19:06:03.97 K4V  6.71£0.16 0.66+0.07 4350-150 0.10-0.02 «A 38 
44 Boo Bb 15:03:46.74 +47:39:15.89 K4V 12.8+0.2 0.55 5035 0.24 46,47,48 
GJ 667 B WDS J17190-3459 B K4V 31 
70 Oph B LHS 459 18:05:27.82 +02:29:39.12 K5V 5.15+0.03 0.73+0.01 4390+200 0.15+0.02 6.2+1.0 1,22 
61 Cygni A HIP 104214 21:06:59.64--38:45:49.90 K5V 3.49740.001_ — 0.67*0 05 43047753 0.15+0.01 6 1,2,26 
61 Cygni B HIP 104217 21:07:00.88 +38:45:21.86 K7V 3.497+0.001 — 0.61*0 05  3949*122 ^ 0.104*0 002 6 1,2,26 
GL 412A HIP 54211 11:05:22.09 +43:31:51.41 M1V  — 4.83-0.03 0.38 3688 0.02070 090 3 2,15,28 
GL 570 B (HIP 73182) 14:57:27.70 -21:25:07.94 M1.5V  5.8640.06 0.5540.05 3345 0.13-£0.09 1,2,37 
GJ 832 LHS 3685, HD 204961 21:33:33.90 -49:00:45.47 M1.5V  4.964--0.001 0.44+0.02 3601 0.03*0 05 9.24 1,2,4 
GJ 667 C 17:19:00.30 -34:59:52.06 M1.5V  7.244+0.005 0.33+0.02 3350450 0.0137+0.0009 >2 31 
Kapteyn’s Star GJ 191 05:11:50.38 -45:02:37.73 M2V  3.9331-:0.0004 0.28+0.02 3742+157 0.016+0.004 115552 1,2,16 
Wolf 1055 HD 180617, GL 752 A 19:16:55.26 +05:10:08.04 M2.5V — 5.9123-0.02 0.4840.02 3534451 0.032-:0.001 40 
AD Leo GJ 388 10:19:35.70 19:52:11.30 | M3V  4.964+0.001 0.42340.012 3477423 0.0236+0.0001 0.40.1 39 
GL 570 C HIP 73182 14:57:27.70 -21:25:07.94 M3V —— 5.86-0.06 0.39+0.03 0.018 37 
Wolf 1061 GJ 628, LHS 419 — 16:30:17.96 -12:40:04.26 M3.5V  4.306+0.001  0.3040.02 33094157 0.01149 005 1,2 
Luyten's Star GJ 273 07:27:25.11+05:12:33.78 M3.5V  5.92-E0.02  0.29230.02 3382+49 0.008749 0065 1,2,43 
55 Cnc B 08:52:40.28 +28:18:54.91 M4.5V 12.48+0.02 0.26+0.02 31874157 0.008+0.002 1,35 
GL 783B HD 191408B, LHS 487 20:11:12.80 -36:06:32.17 M4.5V 6.05+0.03 0.24 0.0009 7.05+0.65 25 
GJ 777B LHS 3509 20:03:37.57 +29:53:45.49 M4.5V — 15.97--0.01 0.21+0.02 3099+157 0.0047+0.0012 1,34 
40 Eri C LHS 25 04:15:19.12 -07:40:15.25 M4.5V 5.04+0.01 0.204+0.006 3100 0.008 6.9+4.7  1,14,23,24 
v And B 01:36:50.16 41:23:25.98 M4.5V 13.47+0.02 0.19+0.02 31594157 0.004479 0012 29 
Proxima Centauri GJ 551 14:29:34.16 -62:40:33.87 M5.5V 1.3012+0.0003 0.1340.02 2900+100 0.001679 0006 4.85 1,17,18 
GL 412B WX UMa, LHS 39 11:05:24.50+43:31:33.27 M6.6V 4.83+0.03 0.095 2700 0.00095 3 27,28 
VB 10 GL 752 B 19:16:57.61+05:09:01.59 M8V  5.912+0.02 0.0881*0 0026 2508183 0.000499+0.000004 1 41,42 
References—1: Stassun et al. (2019), 2: ExoCat-1, Turnbull (2015), 3: Gray et al. (2006), 4: Sanz-Forcada et al. (2010), 5: Huang et al. (2018), 
6: Mamajek & Hillenbrand (2008), 7: Benedict et al. (2022), 8: Gaia Collaboration (2020), 9: Takeda et al. (2007), 10: Gray et al. (2003), 11: 
Delrez et al. (2021), 12: Udry et al. (2019), 13: Pepe et al. (2011), 14: Ma et al. (2018), 15: Mann et al. (2015), 16: Anglada-Escude et al. 


(2014), 17: Faria et al. (2022), 18: Kervella et al. (2003), 19: Ghezzi et al. (2010), 20: Dorren & Guinan (1994), 21: Giidel et al. (1997), 22: 
Eggenberger et al. (2008), 23: Mason et al. (2017), 24: Johnson & Wright (1983), 25: Turnbull & Tarter (2003), 26: Kervella et al. (2008), 27: 
Casagrande et al. (2008), 28: Mann et al. (2015), 29: Rosenthal et al. (2021), 30: Santos et al. (2013), 31: Anglada-Escudé et al. (2013), 32: 
Bidelman (1985), 33: Feng et al. (2019), 34: Hawley et al. (1996), 35: Alonso-Floriano et al. (2015), 36: Gillon et al. (2017), 37: Mariotti et al. 
(1990), 38: Fernandes et al. (1998), 39: Kossakowski et al. (2022), 40: Burt et al. (2021), 41: Pravdo & Shaklan (2009), 42: Pineda et al. (2021), 
43: Astudillo-Defru et al. (2017), 44: Holmberg et al. (2009), 45: Mosser et al. (2008), 46: Zirm (2011), 47: Latkovié et al. (2021), 48: Zasche 


et al. (2009), 49: Philipot et al. (2023), 50: Nordstróm et al. (2004) 


*Significant disagreements in literature are discussed in Appendix A. 


3.2. Chandra 


We retrieved Chandra ACIS imaging observations 
from the Chandra archive’ and reprocessed all evti 
data using the CIAO v4.15 (Fruscione et al. 2006) task 
chandra repro and standard reduction procedures. We 
ran the point source detection algorithm wavdetect on 
each individual exposure to generate a list of X-ray 
source positions and error ellipses. The major and mi- 
nor axes of these error ellipses were increased by a fac- 
tor of five, and all X-ray sources were masked so that a 
background light curve could be extracted for each ob- 
servation. We inspected these background light curves 
for background flaring events. We created good time 


T https:/ /cda.harvard.edu/chaser/ 


intervals (GTIs) using 1c.clean and filtered all event 

data on these GTIs. Finally, we restricted the energy 

range of our resulting “clean” evt2 files to 0.5-7 keV. 
We re-ran wavdetect on the cleaned image to identify 


significant X-ray sources and obtain first-pass estimates 
of their source properties. Due to its exquisite spatial 
resolution (0.496" pixel size) and a typical on-axis PSF 
size less than 1" for the ACIS instruments, Chandra is 
the instrument of choice for resolving crowded fields. 
For stellar systems with multiple components and an- 
gular separations of less than a few arcseconds, Chan- 
dra is currently the only telescope that can be used to 
definitively associate X-ray emission with a given stellar 
component. Figure 2 (right panel) illustrates the spa- 
tial resolving power of Chandra with a 5 ks exposure of 
the triple star system 40 Eri, where the brightest X-ray 


Figure 2. Comparison of the spatial resolution of XMM-Newton (left) and Chandra (right); a 1’ scale bar is shown in both 
images. Left: The GL 412 binary star system with the locations of GL 412A (M1V) and GL 412B/WX Ma (M6.6V) shown by 
cyan circles. The circles have radii of 10" and are separated by —30". A third, unidentified X-ray source (dashed yellow circle) 
is located ~37” from GL 412B. Right: The 40 Eri triple star system with the locations of 40 Eri A (KOV), 40 Eri B (DA white 
dwarf), and 40 Eri C (M4.5V) circled in white. The separation between the B and C components is ^7". X-ray sources are 
found coincident with 40 Eri A and C, while 40 Eri B is undetected. 


Table 3. Summary of X-ray Observations 


XMM-Newton Chandra 
Star # Obs. Exp. Time (ks) # Obs. Exp. Time (ks) 
n Crv 0 ee 4 38.5 
€ Oph 0 mE 1 19.8 
v And 1 5.3 4 58.9 
uL Hor 32 211.3 0 
v Phe 1 16.0 0 
y Pav 1 15.9 0 
B Vir 1 41.1 0 
GJ 1095 0 1 96.2 
LHS 237 1 15.6 0 
t Per 0 1 4.9 
Total 151 2235.7 41 815.0 


NoTE—Select entries are shown to illustrate the table form 
and content. The full machine-readable table (57 rows) is 
available online from the journal. 


source is clearly associated with the position of 40 Eri 
C, a fainter X-ray source is found at the location of 40 
Eri A, and 40 Eri B is undetected. 

If the target star was detected in the cleaned evt2 
image with 250 net counts, we extracted source light 
curves (binned to 100 s) to search for evidence of X- 
ray variability. For sources with 2500 net counts, the 
CIAO tool specextract was used to extract spectra and 
generate RMFs and ARFs. All spectra were binned to 
contain at least 10 counts per bin. 


3.3. Swift 


The X-Ray Telescope (XRT) onboard the Neil Gehrels 
Swift Observatory additionally observed some nearby 
stars. Swiftis optimized to quickly slew to observe tran- 
sient events, and typically obtains short “snapshot” X- 
ray observations (^2 ks or less) with XRT of relatively 
bright X-ray sources. These observations are useful for 
providing flux estimates (or flux upper limits) for stars 
of interest, but often cannot provide the longer-duration 
light curves or higher signal-to-noise spectra that can be 
obtained from XMM-Neuwton or Chandra. We therefore 
use the High Energy Astrophysics Science Archive Re- 
search Center (HEASARC)? to determine which stars 
were observed by Swift but not XMM-Newton or Chan- 
dra, and we use the Swift data products generator? to 
create stacked images of these additional stars and mea- 
sure count rates (or count rate upper limits). The 19 
stars that were observed only by Swift/XRT, but not by 
XMM-Newton or Chandra, are not considered part of 
our primary stellar sample, but we still provide updated 
provide X-ray luminosity estimates in Appendix B. 


4. ANALYSIS PROCEDURE 


We first inspected all X-ray observations for evidence 
of X-ray emission at the location of the target stars. 
For stars that were detected, we first extract light 
curves (binned to 100 s) to search for rapid variabil- 
ity within the observation. For brighter stars (i.e., those 


8 https:/ /heasarc.gsfc.nasa.gov/cgi-bin/W3Browse/swift.pl 
9 https:/ /www.swift.ac.uk/user.objects/index.php 


8 


with 22000 net counts with XMM-Newton or 2500 net 
counts with Chandra), changes in observed count rate 
are readily apparent by visual inspection; for fainter 
stars, a more quantitative assessment of variability is 
needed in addition to visual inspection. 


4.1. Assessing Count Rate Variability 


We use the Anderson-Darling test, which is a non- 
parametric statistical test that is more sensitive to 
short duration variations than the more commonly used 
Kolmogorov-Smirnov (KS) test (Feigelson et al. 2022). 
The Anderson-Darling statistic, A?, is computed as 
(Rahman & Chakrobartty 2004): 


N 
A? -2-N- Y D DIn(^)-& QN 1-208010] 
(1) 


where N is the number of photons observed and t; is the 
arrival time of the i-th bin, scaled to the interval (0,1). 
We use scipy’s anderson routine to compute A? and 
the critical values for the 5% and 1% significance levels 
when compared to a normal distribution. An A? value 
that is larger than a critical value for a given significance 
level is evidence of variability; the typical 596 critical 
value is ~0.75 and the typical 1% critical value is ~1.05 
for the observations considered here. We additionally 
compare the observed light curve to a constant count 
rate and compute a standard reduced x? (x? divided by 
the degrees of freedom) statistic. 

Variability in observed count rate is, to first order, 
driven by a change in X-ray flux incident on the detector. 
In our case, this change in flux is due to changes in the 
intrinsic X-ray luminosity of the target star. However, 
changes in the underlying spectra of the star can addi- 
tionally impart more subtle changes on the light curve, 
as the sensitivity of XMM-Newton and Chandra to X-ray 
photons is energy-dependent. We therefore “flag” seg- 
ments of each light curve according to the type of vari- 
ability observed, hereafter referred to as the VarFlag, so 
that the spectral properties can be measured and an ac- 
curate count rate-to-luminosity conversion factor can be 
established as a function of time during the observation. 
The VarFlags are defined as follows: 


e Q: quiescent periods, when the count rate is low 
and constant within the uncertainties 


e F: strong flaring events, where the count rate is 
above ~80% the quiescent maximum count rate 


e D: decaying periods following flares or descending 
count rates that follow highly elevated count rates; 


these periods are roughly defined as immediately 
following an obvious flare where the count rate is 
below 8096 the maximum flare count rate but still 
above ~150% of the quiescent count rate 


e E: elevated count rate periods, when count rates 
are approximately 5096 higher than the quiescent 
count rate but not obviously associated with a 
strong flaring event 


e R: rising count rate periods that precede a strong 
flare 


Dividing a single observation into a series of sub- 
exposures is a subjective process: we wish to define sub- 
exposures that are short enough to capture unique varia- 
tions seen in the light curve (such as rapid flaring events, 
some lasting only 200-300 s) but long enough so that a 
sufficient number of X-ray counts are available for ade- 
quate spectral modeling (the ~2000 and ~500 count lim- 
its for XMM-Newton and Chandra, respectively). Our 
main goal is to accurately measure the X-ray luminosity 
of the star as a function of time. 

Figure 3 shows an example light curve of Prox- 
ima Centauri (XMM-Newton observation number 
0801880501) and a cumulative distribution function 
(CDF) of photon arrival times compared to a constant 
count rate. The A? statistic, critical values, and x? are 
shown, and the light curve has been color-coded accord- 
ing to the VarFlag assigned to each segment of the light 
curve. Table 4 summarizes the variability metrics for 
every X-ray observation of every star available in the 
archives (the full machine-readable table is available on- 
line from the journal). While some stars in our sample 
have been the targets of extensive observing campaigns 
by XMM-Newton and/or Chandra for many years, ob- 
servational data is quite sparse for other stars. For stars 
with many observations, long-term periodic X-ray vari- 
ability has been detected and is likely the result of coro- 
nal activity cycles (e.g., ; Hor, the a Cen system; Sanz- 
Forcada et al. 2019; Ayres 2023). Similar long-term X- 
ray variability cycles may exist for sparsely observed 
stars in our sample, but will not be detectable in sin- 
gle snapshot observations. Additional observations are 
required to constrain the long-term coronal activity cy- 
cle and/or flaring activity for many stars in the present 
sample. 


4.2. Spectral Modeling 


All X-ray spectral modeling was performed with 
XSPEC v12.11.1 (Arnaud 1996) using standard x? 
statistics. Depending on the number of counts and the 
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Figure 3. Left: Photon arrival time CDF of Proxima Centauri (XMM-Newton observation ID 0801880501, black) compared to 
a constant count rate (dark red, dashed). The A? statistic and critical values are shown in the upper left corner. Right: The 
light curve data (black circles) compared to a constant count rate (dark red, dashed), with sections of the light curve color-coded 
by VarFlag. Quiescent periods are shown in dark blue-violet (*Q"), flaring periods in pink “F”), and descending count rates in 


purple (*D"). 
'Table 4. X-ray Variability Metrics 


Exp. Net 

Star Observation ID Time (ks) | Counts x A? 

u Hor XMM/0673610201 5.3 4230+ 70 1.5 0.49 

B Vir XMM/0044740201 41.1 18450+4160 3.1 0.73 

v And A XMM/0722030101 5.3 940+40 2.1 0.36 
44 Boo XMM/0100650101 18.8 287110+550 3.1 0.64 
GL 620.1A XMM/0822070201 7.8 122204110 1.9 0.41 
k! Cet Chandra/22327 14.0 6600+80 1.7 3.26 
Wolf 1061 Chandra/20163 38.3 25040 0.9 29.78 
Prox Cen XMM/0801880501 20.0 66790+260 35.2 17.55 
GL 412B XMM/0742230101 10.3 320304210 243.1 3.43 
VB10 XMM/0504010101 24.2 660+110 1.1 11.01 


NoTEÉ—Select entries are shown to illustrate the table form and con- 
tent. The full machine-readable table (200 rows) is available online 
from the journal. Net counts are reported for 0.2-15 keV for XMM- 
Newton observations and 0.5-7 keV for Chandra observations. 


overall complexity of the spectra, we used one-, two-, or 
three-temperature APEC models to parameterize the X- 
ray spectra. Since all of the stars are located very close 
to the Sun (within —20 pc), the intervening hydrogen 
column density is expected to be low; we adopt a fixed 
absorbing column Ny (using the TBABS model) of 10!? 
cm"? for all stars (a similar assumption was made by 
Brown et al. 2023). We assume abundances from Wilms 
et al. (2000) for both the plasma emission and inter- 
stellar absorption models. The plasma temperature and 
normalizations were left as free parameters during the 


fit. The normalization of each APEC component is de- 
fined as (neglecting cosmological terms): 


107+ 


where d is the distance to the source, ne and ng are 
the electron and hydrogen column densities (in units 
of cm ?), respectively, and dV is the volume element 
(in units of cm?). The integral f nengdV is the volume 
emission measure (VEM). Once the best-fit model is ob- 
tained, we use the CFLUX convolution model to measure 
the flux on each plasma component. 

While the spectra of some stars required only one or 
two APEC components to achieve a statistically accept- 
able fit, most spectra required a third APEC component. 
We initialized all spectral models with the default (so- 
lar) abundance. For some stars (e.g., 6 Hor, k! Cet) 
the best fit model was insensitive to the exact choice 
of abundance; for these stars the abundance was kept 
fixed at the solar value. If a statistically acceptable fit 
was not obtained with solar abundances, we allowed the 
abundance to be a free parameter of the fit, subject to 
the constraint that the abundance be the same for all 
plasma components. 

For stars that exhibited count rate variability, spectra 
were extracted and modeled for each sub-exposure. Ta- 
ble 5 provides full best-fit parameters for all stars with 
enough counts to perform spectral modeling for all ob- 
servations and sub-exposures (the full machine-readable 
table is available online from the journal) For stars 
with multiple observations or multiple sub-exposures 
with the same VarFlag, we present the average best- 
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Figure 4. Best-fit spectra for all 38 XMM-Newton obser- 
vations of ; Hor (black) with the average best-fit spectrum 
superimposed (red). 


fit model for a given VarFlag in Table 6. For uni- 
formity, we refer to the coolest component of a three- 
component model as “APEC #1,” the hottest compo- 
nent as “APEC #3,” and the intermediate-temperature 
component as “APEC #2.” The two-temperature ther- 
mal plasma models have their individual components 
matched to the three-temperature component with the 
most similar temperature. Figure 4 shows an example of 
the averaged spectrum for ; Hor, which did not exhibit 
significant count rate variability in any of the publicly 
available 38 XMM-Newton observations and was there- 
fore assigned a VarFlag of Q (quiescent). The best-fit 
models for the individual exposures are shown in black, 
with the average best-fit spectrum superimposed in red. 


Figure 5 shows the VEM for each star (calculated from 
the best-fit normalization, defined in Equation 2 above) 
as a function of temperature for the coolest APEC 
component (circles), the intermediate APEC compo- 
nent (squares), and the hottest APEC component (dia- 
monds) using the measurements summarized in Table 6. 
The points are color-coded by stellar spectral type: F 
and G stars are shown in yellow, K stars in orange, and 
M stars in red. In general, the temperatures and VEMs 
of M stars exhibit the largest spread in coronal temper- 
atures relative to more massive stars. The coronal tem- 
peratures of FG and K stars are similar for the coolest 
and hottest APEC components, although VEMs of FG 
stars are roughly an order of magnitude higher than 
those of K stars at similar temperatures. The VEMs 
do not significantly change within a given spectral type 
group, but the intermediate-temperature APEC compo- 
nent of FG stars is on average hotter (by ~2 MK) than 
the intermediate APEC component of K stars. 


The three-component thermal plasma model that gen- 
erally describes the X-ray emission from F and G stars 
is described by plasma temperatures of ~1.3 MK, 5.1 
MK, and 9.4 MK, with corresponding (log) VEMs of 
55.2 [cm?], 55.4 [cm?], and 55.0 [cm], respectively. For 
K stars, the temperatures become 1.2 MK, 3.3 MK, 10 
MK, with (log) VEM values of 54.1 [cm?], 53.9 [cm?], 
and 53.6 [cm?], respectively. For these hotter stars, 
periods of elevated count rates or flaring events differ- 
ent from the quiescent periods in that VEMs increased, 
while the plasma temperatures remained mostly con- 
stant. This result is in sharp contrast to the model fits 
for M stars, which showed more significant variations in 
both temperature and VEM during flaring events. 


4.3. X-ray Luminosity 


In addition to insights about the temperatures and 
VEMs of coronal structures associated with our target 
stars, the best-fit X-ray spectral models provide us with 
a count rate-to-luminosity conversion factor that mit- 
igates against instrumental effects from the extracted 
light curves and gives us a clearer view of the luminosity 
evolution of the target star. Furthermore, we can “un- 
fold” our X-ray spectra (i.e., de-convolve the observed 
spectrum and the energy-dependent detector response 
function) to assess the intrinsic energy emission from 
the target star as a function of energy (or wavelength). 

Figure 6 shows the continuum luminosity of each star 
(in units of Zo) as a function of wavelength during qui- 
escent times. All spectra have been corrected for minor 
foreground absorption, and we use the XSPEC com- 
mand fakeit to extend our best-fit model to a wave- 
length range of ~1-200 A. For nearly all stars, flar- 
ing events are driven by (sometimes dramatic) enhance- 
ments in X-ray luminosity at wavelengths S10 A (en- 
ergies 21.2 keV). Figure 7 shows quiescent and flaring 
spectra for three representative stars (k! Cet, GJ 570A, 
and Prox Cen) from 1-1000 A (top panel), as well as the 
ratio of the flaring and quiescent spectra (bottom panel). 
For both x! Cet (G5V) and GJ 570A (K4V), the flux 
at ~1 A (~12 keV) during flaring times is a factor of 
7-30 larger than during quiescent times, but the fluxes 
at longer wavelengths (2:100 A; energies 0.12 keV) are 
similar. In the case of Prox Cen (M5.5V), however, the 
flaring spectrum is brighter than the quiescent spectrum 
by at least an order of magnitude at all wavelengths, 
with short-wavelength (€ 3 A) emission enhanced by 
more than a factor of 10°. This finding is unlikely to be 
the result of the reduced soft energy responses of Chan- 
dra or XMM-Newton, as the loss of effective area due 


Table 5. All Best-Fit Parameters from X-ray Spectral Modeling 


APEC #1 APEC #2 APEC #3 
Star ObsID VarFlag kT, Ni kT> No kT3 N3 logLx abund (x?/dof) 
(keV) (1074) (keV) (1074) (keV) (1074) [erg s~1] 

n Crv 14474 Q 0.164997 sat? 0.50+0.08 4.5 20 28.433:-0.060.32*026 — 1.14 
B Vir 0044740201 Q  0.0020.01 4.7*L2 0.3740.01 2.4+0.1 28.08+0.02 1 (fixed) ^ 2.48 
x! Ori 0111500101/0-lOks E 0.10+0.02 39.71387 0.49-:0.0352.9* 10,5 0.84* 1-00 12.1*7$ 28.881903 0.311665 — 1.31 
B Com  0148680101/0-bks Q  0.072-0.0223.5* 1145 0.394-0.02 3.4+0.3 28.11+0.02 1 (fixed) 1.55 
44 Boo 0100650101/0-6ks D  0.25+0.02 23.87¢¢ 0.50*015 9.8134 0.917505 275.3*12228.672-0.010.09-0.01 — 1.45 
GL 620.1A ^ 0822070301 Q 011992 is UO 1.60 


—0.01 


9.5433 0.46+0.03 1.4 


56819. 28.82 


MT +0.09 
+£0.03 0.43* 0-09 


NoTE—Select entries are shown to illustrate the table form and content. The full machine-readable table (360 rows) is available 


online from the journal. 


Table 6. Average Best-Fit Parameters from X-ray Spectral Modeling 


APEC #1 APEC #2 APEC #3 
Star VarFlag kT Ni kT> N2 kT3 N3 logLx abund (x? /dof) 
(keV) (1074) (keV) (10-4) (keV) (10-4) [erg s-!] 

gow Q 0.204912 Baty) 0.551099 3.1122 28.41+9:08 Qar oe 1.31 
vAndA Q 0.29+0.03 0.80.1 27.57+0.03 1 (fixed) 0.91 

t Hor Q  0.09+0.03 7.7+1.1 0.42+0.08 2.9+0.2 0.71+0.14 1.6+0.2 28.75*0.10 1 (fixed) — 1.27 
BVir Q  0.09+0.01 4.7*12 0.37+0.01 2.4+0.1 28.08+0.02 1 (fixed) 2.48 
x! Ori Q — 0.08+0.03 258+228 0.42+0.07 28.6+9.2 0.69+0.06 21.3 0.12 1.55 


+0.07 
8.2 28.94* 007 0.204 


NoTE—Select entries are shown to illustrate the table form and content. The full machine-readable table (76 


rows) is available online from the journal. 


to contamination build-up occurs most significantly at 
energies below 1 keV (e.g., Grant et al. 2024). Never- 
theless, we use detector response functions specific to 
the orbital cycle in which the individual observations 
were obtained to de-convolve the detector response as 
accurately as possible. Plots of individual stellar X-ray 
spectra, color-coded by VarFlag, are discussed in Ap- 
pendix A. 

A question that our analysis can elucidate is: For what 
fraction of time is a star observed to be emitting at a 
given X-ray luminosity? While this is related to the 
flaring frequency and the typical peak luminosity and 
duration of flares, our goal in this work is not to directly 
characterize flaring events or the physical mechanisms 
that cause them. Rather, we wish to predict the X-ray 
environments in which HZ planets reside. 

'To do this, we use our best fit spectral models to first 
convert the observed count rates to luminosities. For 
stars that are too faint to have their spectra directly 
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modeled, we use the spectral model of the star with the 
closest spectral type and physical parameters to convert 
the observed count rates into luminosities. For stars 
that are detected but too faint for reliable light curve ex- 
traction, we simply estimate their time-averaged X-ray 
luminosity; for un-detected stars we calculate 3e lumi- 
nosity upper limits. We then apply Astropy’s (Astropy 
Collaboration et al. 2013, 2018) bayesian_block algo- 
rithm to the light curves; the “fitness” parameter is used 
to define the measurement uncertainties. The Bayesian 
block algorithm (Scargle 1998; Scargle et al. 2013) is a 
segmentation technique, which aims to optimally split 
(in this case) time series data into “blocks” such that 
each block is statistically different from its neighboring 
block. An example Bayesian blocked light curve for «t 
Cet is shown in Figure 8. 

The output of the bayesian_block routine is a list of 
time bins and the luminosity the star was emitting at 
during each bin. From these lists, we measure the frac- 
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Figure 5. Best fit VEMs and temperatures for thermal plasma components from XSPEC modeling. VEMs and temperatures 
from APEC #1 are shown in circles, those from APEC #2 are shown in squares, and those from APEC # are shown in 
diamonds. Values measured for F and G stars are shown in yellow, K stars in orange, and M stars in red. M stars exhibit the 
largest variations in both VEM and temperature across all three APEC components. 


n Crv (F2V) —— 55 Cnc A (KOIV/V) 
— v And (F8V) —— 70 Oph (KOV--K5V) 
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Figure 6. Average continuum quiescent spectra for each target star (0.5-150 A, corresponding to an energy range of ~25-0.08 
keV). Spectra have been corrected for minor foreground absorption. Stars are color-coded by spectral type, with earlier-type 


stars shown in yellow and later-type stars shown in dark red. 


tion of the time that each star was observed to be fainter 
than a given Lx or Lx / Lys; value. The results are shown 
in Figure 9. For comparison, Lx of the modern-day Sun 
is 6.31x107° erg s^! (logLx/Lyoy = -6.78) when quiet 
and 8.04x 10?" erg s~! (logLx / Lyoy = -5.68) when active 
(Linsky et al. 2020), while ~4 Gyr ago the Ly of the 
Sun was ~75% the modern value and its Lx was a fac- 
tor of ~2-3 higher (Obridko et al. 2020). On this plot, 
a perfectly non-variable star would be a straight verti- 


cal line, while stars that experience significant changes 
in X-ray luminosity (i.e., Proxima Centauri) will show 
a pronounced horizontal extent. We measure the lumi- 
nosities at which the stars spend less than 25%, 50%, 
70%, and 90% of the observed time (which we refer to 
as Lo5, Lso, L75, and Lgo, respectively). These values 
are summarized in Table 7; we also include the time- 
averaged luminosities (with uncertainties) of faintly de- 


—— kl Cet (G5V) 
—— GJ 570 A (KAV) 
—— Prox Cen (M5.5V) 


109 101 10? 103 
wavelength (A) 


Figure 7. Top: Quiescent (solid lines) and flaring (dashed 
lines) continuum spectra of three representative stars from 
our sample (&! Cet in yellow, GJ 570 A in orange, and Prox 
Cen in red). Bottom: The ratio of flaring spectrum to quies- 
cent spectrum, showing a pronounced enhancement in X-ray 
emission at «10 À during flaring periods for all three stars. 
The flaring spectrum of Prox Cen remains an order of mag- 
nitude brighter than the quiescent spectrum even at longer 
wavelengths. 
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Figure 8. The light curve of k Cet from XMM-Newton ob- 
servation 0111410101. The MOSI count rates were converted 
to luminosities (gray circles) using the best-fit spectral mod- 
els from three sub-exposures (from 0-14 ks, 14-25 ks, and 
25-39 ks). The Bayesian blocked light curve is superimposed 
in blue. 
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Figure 9. Top: The fraction of time that a star is observed 
to be fainter than a specific luminosity. Bottom: The frac- 
tion of time that a star is observed to be below a Lx/Lpo1 
value. Stars are color-coded by spectral type, with earlier- 
type stars shown in yellow and later-type stars shown in dark 
red. The 25th, 50th, 75th, and 90th percentile luminosities 
and Lx/L»o ratios are indicated by blue dashed lines (val- 
ues for each star summarized in Table 7). The shaded gray 
region indicates typical ranges for the current Sun, and the 
blue shaded regions indicate approximate ranges for the early 
Sun (~4 Gyr ago; see text). 


tected stars and the 30 upper limits for non-detected 
stars in the L5o column. 


5. DISCUSSION 


We compared the 25th and 90th percentile Lx / Lsoi 
ratios for stars in our sample to the current minimum 
and maximum ratios for the Sun (Linsky et al. 2020), 
and find that eight stars that were detected in the X-ray 
observations at high significance exhibit Lx/Lypo ratios 
similar to that of the modern Sun: £ Oph (F2V), 8 Vir 
(F9V), 47 UMa (GOV), 8 Com (GOV), 55 Cnc A (G8V), 
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Table 7. X-ray Luminosities and Lx/Lpo1 Ratios (or 30 Upper Limits) for All Sample Stars 


logLos  logLso  logLzs  logLeo log(L25/Lpo1) log(Ls50/Lpo1) log(L75/Lo1) log(Lo0/Lpo1) 
Star [erg s71] [erg s] [erg s+] [erg s- 1] 

n Crv 28.19 28.39 28.47 28.59 -6.06 -5.86 -5.78 -5.66 
€ Oph 27.47 27.64 27.77 28.03 -6.73 -6.56 -6.43 -6.17 
v And A 27.27 27.55 27.12 27.87 -6.84 -6.56 -6.39 -6.24 
B Vir 27.90 27.99 28.05 28.15 -6.26 -6.17 -6.11 -6.01 
B Hyi 26.23 26.53 26.94 27.37 -7.92 -7.62 -7.21 -6.78 
k! Cet 28.63 28.71 28.78 28.88 -4.89 -4.81 -4.74 -4.67 
HD 140901 27.62 27.92 28.05 28.23 -6.05 -5.75 -5.62 -5.44 
ô Eri 26.39 26.53 26.69 26.89 -7.72 -7.58 -7.42 -7.22 
GL 570 A 27.15 27.36 27.52 27.68 -5.78 -5.57 -5.41 -5.25 
61 Cyg A 26.75 26.92 27.08 27.28 -6.00 -5.84 -5.68 -5.56 
61 Cyg B 26.69 26.87 27.04 27.42 -5.89 -5.71 -5.54 -5.35 
AD Leo 28.42 28.51 28.69 29.36 -3.54 -3.45 -3.27 -2.98 
Proxima Centauri 26.70 26.93 27.71 28.63 -4.09 -3.86 -3.08 -2.33 


NoTE—Select entries are shown to illustrate the table form and content. The full machine-readable table (63 rows) is available online from 


the journal. 


GL 892 (K3V), GL 783 (K3V), and GL 183 (K3V). 
Three of these stars (8 Vir, 8 Com, and GL 183) are 
younger than the Sun, while four are older (47 UMa, 55 
Cnc A, GL 892, and GL 783) and one (£ Oph) does not 
have an age estimate available in the literature. Twelve 
stars have Lx /Lyo) ratios consistent with those of the 
early Sun: 7 Crv (F2V), 8 Com, GL 311 (G1.5V), c! 
Ret (G2V), GL 327 (G3V), HD 140901 (G7IV), the un- 
resolved binary 70 Oph (KOV+K5V), LHS 1875 (K3V), 
GL 570A (K4V), 61 Cyg A (K5V) and B (K7V), and 
the unresolved GL 570 BC pair (M1.5V--M3V). We note 
that one star, 8 Com, exhibits a range in Ly / Ly that 
fits within the narrow overlap region between the mod- 
ern and early Sun (logLx/Lyaz -6.4 to -5.7). An ad- 
ditional 11 stars were not detected or only marginally 
detected and have Lx/Lyo upper limits lower than the 
modern day minimum solar Lx/Lyo: v Phe, LHS 237, 
LHS 208, & Per, p CrB, GL 672, 51 Peg, GJ 777A, 
82 Eri, GL 451A, and Kapteyn's Star. For roughly a 
dozen stars in our sample, the observed Lx (or Lx up- 
per limit) translates into X-ray surface fluxes below the 
^-10* erg s~! cm7? minimum observed in the quiet Sun 
and other “coronal hole" stars (Schmitt 1997). A closer 
look at these very X-ray faint stars will be the subject 
of an upcoming study (Binder et al., in preparation). 
In Table 8, we summarize these X-ray-detected Solar- 
like stars (their spectral types and approximate ages are 
taken from Table 2) and the relevant solar Lx / 1.) com- 
parison (modern or early Sun). We additionally pro- 
vide the “tier” each star was assigned in Mamajek & 
Stapelfeldt (2024), indicating the priority which each 
star should be considered for an HWO direct imaging 


survey (with A indicating the highest priority group, 
and tiers B and C requiring further study). 

Any HZ planets that may exist around the stars listed 
in Table 8 are currently experiencing high energy radi- 
ation environments that are similar to or less hostile 
than the modern Earth. Out of the 229 unique stars 
we identified as potential targets for future direct imag- 
ing surveys, only ~25% had ever been imaged with ei- 
ther XMM-Newton or Chandra (this fraction increases 
to ~30% when Swift imaging is included), and out of the 
57 stellar systems with high quality X-ray observation 
we find that 28 systems (~50%) have Lx/Lwo ratios 
similar to or less than that of either the modern or early 
Sun. These stars span the full range of spectral types 
and ages considered here, and represent all three tiers 
of the Mamajek & Stapelfeldt (2024) ranking system. 
It is therefore likely that there are an additional ~100 
potential direct imaging target stars with HZs with cur- 
rent X-ray radiation environments similar to or even less 
hostile than that of the Sun. We note that threshold 
Lx /Lpoiratios below that of the modern or early Sun 
could be a conservative condition for finding habitable 
planets. Although, the past evolution of the star would 
also have to be taken into account. Detailed coupled 
geological and atmospheric evolution analyses would be 
most informative for predicting the targets with the best 
chances of maintaining potentially habitable conditions 
(e.g., Krissansen-Totton & Fortney 2022). A systematic 
X-ray survey of these as-of-yet unobserved nearby stars 
is needed to inform the selection of the best targets for 
direct imaging surveys. 
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Table 8. Potential HWO Target Stars with Solar-Like Lx / Lsoi Ratios 


Star Name Spectral Type Age (Gyr) Sun Comparison HWO Tier 
n Crv F2V <2 early C 
€ Oph F2V unknown modern C 
v Phe F9V ^3 «modern min B 
B Vir F9V ~3 modern C 
LHS 237 F9V unknown <modern min none 
u Per F9.5V ~4 <modern min A 
LHS 208 GOV ~3 <modern min B 
p CrB GOV ~11 <modern min B 
GL 672 GOV ~12 <modern min C 
47 UMa GOV ~6.5 modern A 
B Com GOV <2 modern,early A 
GL 311 G1.5V <2 early B 
C! Ret G2V «2 early A 
51 Peg G2V ~T <modern min none 
GL 327 G3V <2 early C 
GJ 777A G7IV/V ~13 <modern min B 
HD 140901 GTIV/V ^3 early C 
82 Eri G8V ^6 «modern min B 
70 Oph AB KOV+K5V ~6 early B 
55 Cnc A KOIV/V ~9.5 modern C 
GL 451A KIIV ^5 «modern min C 
GL 783 K2.5V ~T modern B 
GL 892 K3V ~11 modern A 
GL 183 K3V ~2 modern B 
LHS 1875 K3.5V <3 early C 
GL 570A K4V <3 early A 
61 Cyg AB K5V+K7V ~6 early A 
GL 570 BC M1.5V+M3V <3 early none 
Kapteyn’s Star M2V ~11 <modern min none 


In general, M stars exhibit significantly more X-ray 
variability than earlier type stars and spend a signif- 
icant amount of time at Ly/Lyo) ratios dramatically 
larger than that of the Sun. We can supplement our 
sample of Sun-like and low mass stars with X-ray ob- 
servations of an 23 additional M- and K-type stars from 
the MUSCLES/Mega-MUSCLES survey (Brown et al. 
2023). Brown et al. (2023) found a major decrease in 
the X-ray-to-bolometric luminosity ratio (Lx /Lo1) with 
increasing effective temperature. We compute Lx / Lyai 
(or upper limits, in the case of X-ray non-detections) 
for all stars in our sample using Ds) from Table 7 
and Lpo) from Table 2. In Figure 10, we plot the 
(logarithmic) Lx/Lyo ratio as a function of Teg (also 
from Table 2) for our full sample (colored points), with 
the MUSCLES/Mega-MUSCLES sample shown in gray. 
Our sample extends this inverse relationship between 
Lx/Lya and Teg by ~1000 K in Teff, with Lx/Lya ra- 


tios that are about an order of magnitude lower than 
observed for K and M stars. 

There is a small clustering of stars with Tog ~ 5200 — 
6200 K with significantly higher Lx /Lpo ratios than the 
majority of comparable stars, likely due to young stel- 
lar ages. Chromospheric activity in solar-type stars is 
observed to decrease rapidly after ~2-3 Gyr (Pace & 
Pasquini 2004; Zhang et al. 2019), after which activity 
levels remain mostly constant. We indicate the approxi- 
mate age of each star in Figure 10: young stars ($2 Gyr) 
are shown in red, older stars (Z2 Gyr) are shown in pur- 
ple, and stars without age estimates or ambiguous ages 
are shown in dark blue. Most of these high-Lx / Loi; 
high-T. stars are indeed believed to be consistent with 
the younger ages where each spectral type emits elevated 
high energy emission (Richey-Yowell et al. 2023). There 
are three stars with ambiguous ages that lie within this 
region of the diagram (x! Ori, 55 Cnc A, and GL 451A); 
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Figure 10. The Lx/Lpo ratio as a function of Tig. Stars 
in our sample with firm X-ray detections are shown in cir- 
cles, while upper limits are shown for stars that were not 
X-ray detected (downward pointing triangles). Stars that 
are believed to be young (<2 Gyr) are shown in red, older 
stars (22 Gyr) are shown in purple, and stars without age 
estimates or with ambiguous/discrepant age estimates are 
shown in dark blue. The MUSCLES/Mega-MUSCLES sam- 
ple (Brown et al. 2023) is shown in gray. The yellow circles 
and shaded area indicate the quiescent and flaring Lx /Lpol 
ratios for the Sun (Linsky et al. 2020). 


the observed X-ray emission of these stars suggests that 
younger age estimates may be preferred for these stars. 

In Figure 11 we again show the Lx/Lpo ratio as a 
function of Tep for stars in our sample, but we use 
Los, Lso, Lrs and Loo to compute the Lx /Lpol ratio. 
This allows us to assess the degree to which intrinsic 
X-ray variability contributes to the scatter in Lx /Lpol 
observed in Figure 10. The ratio of L75/ L55 provides an 
estimate of how much a star varies “normal” about the 
median Lx. We find that the variability of the youngest 
stars in our sample exhibit the smallest dynamic range, 
with L75/L25 ~ 1.3 for stars with Teg 25000 K (for the 
coolest young star in the plot, AD Leo, Lz5/ Las ~ 1.8). 
The Lz5/ Ls ratio for older stars shows more dispersion, 
varying by up to a factor of ~3.5. 
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Figure 11. Top: the Lx/Lpa ratio as a function of Teg as- 
suming Le5 (black), Lso (blue), Lz5 (green), and Lgo (dark 
yellow). The yellow circles and shaded area indicate the qui- 
escent and flaring Lx /LZyo1 ratios for the Sun (Linsky et al. 
2020). Bottom: the ratio of L75/L25 as a function of Ter, 
color-coded by age (as in Figure 10: young stars shown in 
red, old stars shown in purple, and stars with ambiguous 
ages shown in dark blue). 


6. CONCLUSIONS 


We have analyzed —3 Msec of high quality XMM- 
Newton and Chandra observations of 57 nearby stellar 
systems, as well as snapshot Swift images of an addi- 
tional 19 stellar systems, to provide a uniform catalog 
of the X-ray properties of potential target stars for fu- 
ture direct imaging surveys. All X-ray light curves and 
spectroscopic data is publicly available on Zenodo (DOI: 
10.5281/zenodo.11490574). The X-ray environments in 
the HZs around F, G, and early K stars are more likely 
to be similar to those that the Earth has evolved in, and 
we identify 29 stars with current Lx/Lypo ratios similar 
to or less than that of the Sun. Our results are consis- 
tent with Richey-Yowell et al. (2023), who found that 
the X-ray fluxes of late M stars can be a factor of ~3-15 
times larger than for K stars. In addition to high aver- 
age Lx/Lya ratios, M stars exhibit more frequent and 
dramatic flares than earlier-type stars. We note that, of 
the 229 candidate target stars for direct imaging, only 
~30% have been observed by XMM-Newton, Chandra, 
or Swift. 
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APPENDIX 


A. DISCUSSION OF INDIVIDUAL STARS 


In this appendix we provide additional context for 
each stellar system (e.g, when there is a notable dis- 
crepancy in the published literature with regards to a 
given physical parameter) and summarize the known 
or candidate exoplanets associated with each star. We 
then present the results of our X-ray analysis for each 
stellar system. Table 9 summarizes the observations of 
faint/non-detected stars, including if the star was de- 


tected in a given observation or not, the count rate (or 
upper limit on the count rate) and detection significance, 
and the “match star" that was used to convert count 
rates to luminosities. 

Figure 12 shows X-ray images of resolved and unre- 
solved binary star systems included in our study. In- 
dividual images are discussed further for each star be- 
low. Figure Set 13 shows the light curves extracted for 
all stars from all observations in which they were de- 
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tected, with sub-exposures colored by VarFlag. Figure 
Set 14 shows the average best-fit spectra for all suffi- 
ciently bright stars. For stars that exhibit flaring events, 
we also show the ratio of the flaring-to-quiescent spectra 
as a function of wavelength. 


A.l. m Crv 


There are no significant discrepancies in the stellar 
physical parameters of 7 Crv. Although there are no 
known or candidate exoplanets, the system is known to 
host two debris disks: one at a distance of ~160-180 AU, 
and a hotter disk at ~3-7 AU (Wyatt et al. 2005, 2007; 
Smith et al. 2008). 

n Crv was observed four times with Chandra. All ob- 
servations were taken with ACIS-S in 1/4 subarray mode 
for ~ 10 ks. In all observations, 7 Crv was detected 
with ~1000 net counts (0.5-7 keV), but no significant 
variability was observed in the X-ray light curves (see 
Figure Set 13). The spectrum of 7 Crv is well-described 
by a two-component thermal plasma model, with no sig- 
nificant changes in the temperatures or normalizations 
across the different observations. 


A.2. € Oph 


There are no significant discrepancies in the stellar 
physical parameters of £ Oph, and the star does not 
host any no known or candidate exoplanets. £ Oph was 
observed once with Chandra/ACIS-S for 19.8 ks and 
robustly detected with ~160 counts (80 significance). 
'There are hints of low-level variability in the X-ray emis- 
sion of € Oph, but the relatively low number of net 
counts prohibits spectral modeling. 


A.3. v And 


v And A is an F8V dwarf in a presumed binary system 
with a proper motion companion v And B, an M4.5V 
dwarf (also referred to as v And D). The two stars have 
an angular separation of —0.9' (corresponding to ~750 
AU). There are no significant discrepancies in the phys- 
ical parameters of either star. v And A hosts three 
confirmed exoplanets (v And b, c, and d) which were 
discovered via RV measurements. The masses of v And 
b, c, and d are My sini ~0.67 Mjup, ~0.67 Myup, ~2 
Mjup, and ^4 Mjyp, respectively, with orbital periods 
(orbital distances) of ~4.6 days (0.06 AU), ~241 days 
(0.82 AU), and ~1282 days (2.5 AU), respectively (all 
planetary data taken from Rosenthal et al. 2021). 

The v And system was observed four times with 
Chandra/ ACIS-S in 1/8 subarray mode (for ~15 ks per 
observation) and once with XMM-Newton for 5.3 ks. v 
And A was detected in all five observations with suffi- 
cient counts to extract light curves (see Figure Set 13) 


and with enough counts in Chandra observations 10976 
and 10977 to enable spectral modeling. There is no ev- 
idence for strong X-ray variability in any observation, 
and both Chandra spectra were consistent with a single- 
temperature thermal plasma model with a temperature 
of ~0.3 keV (~3.4 MK; see Table 6). v And B was out of 
the field of view of the Chandra observations. However, 
the star was within the field of view of the XMM-Newton 
PN detector but fell just on the gap between two CCDs 
and was not detected (see Figure 12). 


A.4. & Hor 


, Hor is characterized either as a late F-type (F9V; 
Gray et al. 2006) or early G-type (GOV; Turnbull 2015) 
star. There is additional disagreement on the age of 
the system: Turnbull (2015) quotes an age of 2.72 Gyr, 
while Sanz-Forcada et al. (2010) suggest a much younger 
age of ~0.47 Gyr (see also the detailed analysis of Sanz- 
Forcada et al. 2019, where an age of 600 Myr is adopted). 
We adopt a spectral type of F9V and an age of 0.47 
Gyr in this work. ; Hor hosts one confirmed exoplanet 
that was first detected via radial velocity (RV) measure- 
ments by Kürster et al. (2000). It is a giant planet with 
My sini~2.3 Mj, and an orbital period of ~300 days 
(corresponding to an orbital distance of ~0.9 AU; Stas- 
sun et al. 2017). 

t Hor is one of the most thoroughly observed stars 
in our sample, with 32 publicly available XMM-Newton 
observations taken between 2011 and 2018. A detailed 
study of the coronal activity and structure of the star is 
presented in Sanz-Forcada et al. (2019), who combined 
the XMM-Newton monitoring with observations from 
TESS and the STIS instrument on the Hubble Space 
Telescope. 'They corroborate a 1.6 year X-ray activity 
cycle for the star, although . Hor does not exhibit the 
same type of dramatic changes in luminosity as seen in 
other stars in our sample (see Figure 9). Our X-ray 
analysis of . Hor is consistent with Sanz-Forcada et al. 
(2019), although we do not undertake a more extensive 
study of the abundances of the star. Light curves of & 
Hor are shown in Figure Set 13, and the best-fit spec- 
trum is described in Table 6. 


A.5. v Phe 


v Phe is characterized either as a late F-type (F9V; 
Gray et al. 2006) or early G-type (GOV; Turnbull 2015) 
star. There are no confirmed or candidate exoplanets 
in the system, and the star was undetected in a 16 ks 
XMM-Newton observation. 


A.6. y Pav 


y Pav is a known F-type star, with the specific spec- 
tral type ranging from F9V (Gray et al. 2006) to F7V 
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ZA GL.620.1A 
: 
Q : 


WD 1620-391 


Wolf 1055 


(g) (h) (i) 

Figure 12. Images of resolved and unresolved binary systems with XMM-Newton and Chandra. (a) XMM-Newton/PN image 
of v And. The locations of the F8V star (A) and the M4.5V star (B) are shown by cyan circles (both with radii 10"). The stars 
are separated by ~0.9’. (b) XMM-Newton image of GL 620.1A and the white dwarf WD 1620-391. No X-rays are detected at 
the location of WD 1620-391 in either observation. (c) The o Cen binary system as seen by the PN camera on XMM-Newton. 
The positions of the G2V star (A) and K1V star (B) are shown in blue circles (both with radii of 5" and separated by ^18"). 
(d) The € Boo binary system as seen by the MOS2 camera on XMM-Newton. The locations of the G8V star (A) and K4V 
star (B) are shown in cyan circles (both with radii of 5"). (e) The 70 Oph binary star system as seen by the PN camera on 
XMM-Newton. The positions of the K0V star (A) and the K5V star (B) are shown in blue circles (both with radii of 5" and 
separated by —3.7"). (f) The GL 783 binary star system as seen by the PN camera on XMM-Newton. The positions of the K3V 
star (A) and the M4.5V star (B) are shown in cyan circles (both with radii of 5" and separated by ~2.6’. (g) Chandra image 
of the GL 570 system. The primary GL 570A (K4V) is easily resolved from the binary red dwarf pair GL 570BC (M1.5V and 
M3V, respectively). Cyan circles have a radius of 10". The brown dwarf GL 570D is not detected in X-rays. (h) The 61 Cyg 
binary system as seen by the PN camera on XMM-Newton. The positions of the K5V star (A) and K7V star (B) are shown in 
cyan circles (both with a radius of 10" and separated by —12"). (i) Chandra image of Wolf 1055 and VB 10 (both circles have 
radii of 10" and are separated by ~1.2’). 

Some images have been smoothed for display purposes only; white contours indicate the intensity distribution of X-ray counts 
on the detector. 
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Figure 13. Left: Photon arrival time CDF for a Cen (black) compared to a constant count rate (dark red, dashed). The 
A? statistic and critical values are shown in the upper-left corner. Right: The light curve data (black circles) compared to a 
constant count rate (dark red, dashed). Quiescent periods are shown in blue-violet, flaring periods in pink, elevated count rate 
periods in blue, and descending periods in purple. The complete figure set (193 images) is available in the online journal. 
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Table 9. Summary of Faint Star Observations 


Star Observation ID Detected? Count Rate (ct s-!) Significance (c) Match Star 
€ Oph Chandra/27852 Yes (8.141.0)x 107? 8.0 9 Crv 
v Phe XMM/0206540101 No «4.3x10-? ee L Hor 
+ Pav XMM/0670380101 Yes (1.140.4) x 107? 2.8 t Hor 
GJ 1095 Chandra/4199 Yes (2.140.7)x 1074 2.9 ı Hor 
LHS 237 XMM/0840210501 No <3.7x 107? es ı Hor 
LHS 208 XMM/0865400201 No <1.4x107? GL 620.1A 
LHS 208 Chandra /22293 No <5.7x1074 GL 620.1A 
t Per Chandra/12338 No «1.2x107? ee L Hor 
p CrB Chandra/12396 No «3.8x10-? ee L Hor 
GL 672 Chandra/12397 No «2.1x10-? tee GL 620.1A 
47 UMa XMM/0304203401 Yes (3.90.3 x 1072 12.1 GL 620.1A 
B Hyi Chandra/12337 Yes (5.140.6 x107? 12.5 GL 620.1A 
B Hyi XMM/0006010401 Yes (5.60.4 x107? 14.0 GL 620.1A 
18 Sco XMM /0303660101 Yes (2.20.3 x107? 7.3 GL 311 
18 Sco Chandra/12393 Yes (5.5+1.1)x 1074 5.0 GL 311 
u Ara XMM/0551021001 No <1.6x107? GL 327 
u Ara XMM/0551023101 No «2.2x10-2 exa GL 327 
HD 136352 XMM/0884680201 Yes (4.141.5 x 1073 2.7 Kl Cet 
51Peg = XMM/0551020901 No <2.3x 1078 k! Cet 
51 Peg Chandra/10825 No <3.1x1073 Kl Cet 
GJ 777A XMM/0304201101 No «1.2x107? es ó Pav 
GJ 777A | XMM/0304202601 No «T.8x10-? ee ô Pav 
HD 140901 Chandra/13769 Yes (8.92:0.9) x 1073 11.0 ó Pav 
82 Eri XMM/0670380601 No <1.2x10~? 55 Cnc A 
82 Eri Chandra/22292 No «8.1x1074 55 Cnc A 
GL 451A Chandra/9931 No «4.9x10-4 GL 183 
T Cet XMM/0670380501 Yes (1.20.2) x 107? 7.0 55 Cnc A 
7 Cet Chandra/1886 Yes (2.4+0.3) x 1073 8.0 55 Cnc A 
6 Eri XMM/0205720101 Yes (7.81.1) x 10-3 7.3 55 Cnc A 
40 Eri A Chandra/13644 Yes (1.42-0.2) x 1072 7.7 GL 117 
GL 412 A  XMM/0742230101 Yes (5.740.2)x 10-2 28.5 GJ 832 
GJ 667 C Chandra/17317 Yes (5.72:0.9) x 1073 6.3 GJ 832 
GJ 667 C  Chandra/17318 Yes (3.140.5)x 1074 6.2 GJ 832 
Kapteyn Chandra/merged No «1.2x10-4 Wolf 1055 
Luyten Chandra/20164 Yes (3.540.7) x 107? 5.0 AD Leo 
v And B XMM/0722030101 No «3.2x10-? Ue 40 Eri C 
VB 10 XMM/0504010101 Yes (2.740.5 x 1072 5.4 GL 412B 
VB 10 Chandra/616 Yes (1.20.8 x 1073 1.5 GL 412B 
VB 10 Chandra/7428 Yes (1.60.4) x 103 4.0 GL 412B 
(Turnbull 2015). There are significant discrepancies in A.T. B Vir 


the estimated age of y Pav: Holmberg et al. (2009) quote 
an age of ~1 Gyr based on photometric observations and 
comparison to stellar evolution models, while and aster- 
oseismic analysis by Mosser et al. (2008) yields an age of 
7.25 Gyr. There are no known or candidate exoplanets 
orbiting y Pav. The star was only marginally detected 
(2.80) in a 25.9 ks XMM-Newton observation. The low 
X-ray luminosity of y Pav is consistent with the older 
age estimate of Mosser et al. (2008). 


There are no significant discrepancies in the stellar 
physical parameters for @ Vir, and the star does not 
host any currently known or candidate exoplanets. 8 Vir 
was observed once with XMM-Newton (observation ID 
0044740201) and was detected with ~18,450 net counts 
on the PN camera (thick filter). The X-ray light curve 
does not show evidence of significant variability (see Fig- 
ure Set 13). The spectrum is well described by a two- 
temperature thermal plasma model (see Table 6). 
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Figure 14. Top: The best-fit continuum spectra for œ Cen 
color-coded by VarFlag. Quiescent periods are shown in dark 
blue, elevated count rate periods in light blue, flaring periods 
in pink, periods of descending count rates in purple, and 
periods of rising count rates in red. Bottom: The ratio of the 
flaring-to-quiescent spectrum as a function of wavelength. 
The complete figure set (18 images) is available in the online 
journal. 


A.8. GJ 1095 


'The spectral type of GJ 1095 is given as F9V by Gray 
et al. (2003) or GOV by Turnbull (2015). There are no 
known or candidate exoplanets orbiting GJ 1095. The 
star was only marginally detected (2.90) in a 96.2 ks 
Chandra observation. 


A.9. LHS 287 


There are no significant discrepancies in the stellar 
physical parameters of LHS 237, and no age estimate 
is currently available in the literature. LHS 237 is a 
triple-star system: the inner spectroscopic binary has 
an orbital period of ~10 years (the angular separation 
of ~0.3” is unresolvable by XMM-Newton; Tokovinin 
2012), and the third component is one of the coolest 
known white dwarfs (NLTT 18141 = GJ 288B) located 
7-14.5' from the central binary (Holberg et al. 2013). 
'There are no known or candidate exoplanets in the sys- 
tem. 

LHS 237 was serendipitously imaged by XMM-Newton 
for 15.6 ks in observation 0840210501 (the intended tar- 
get for the observation was a Wolf-Rayet star). The in- 
ner spectroscopic binary is not detected, and the white 
dwarf is outside the XMM-Newton field of view. 


A.10. & Per 


'There are no significant discrepancies in the stellar 
physical parameters for ; Per. The spectral type is either 
F9.5V (Gray et al. 2003) or GOV (Turnbull 2015). There 
are no known or candidate exoplanets in the system. 4 
Per was observed for 4.9 ks with Chandra in ACIS-S 1/4 
subarray mode, but the star was not detected. 


A.ll. B Hyi 


'There are no significant disagreements in the stellar 
physical parameters of 8 Hyi, and the star does not host 
any currently known or candidate exoplanets. 6 Hyi 
was observed once with Chandra/ ACIS-S in 1/4 subar- 
ray mode and once with XMM-Newton. The star was 
detected in both observations with sufficient counts to 
extract light curves, but the star is too faint for accurate 
spectral modeling to be performed. There is no strong 
evidence for X-ray variability in either the Chandra or 
XMM-Newton light curves of 8 Hyi (see Figure Set 13). 


A.12. B Com 


'There are only minor differences in the reported age of 
8 Com: Takeda et al. (2007) report an age «1.12 Gyr, 
while Turnbull (2015) report an age of 1.7 Gyr. All other 
spectral parameters are taken from Stassun et al. (2017). 
'There are currently no known or candidate exoplanets 
orbiting 8 Com. 

8 Com was robustly detected (with ~38,000 net 
counts) in a 57.2 ks observation by XMM-Newton. 
We extracted a light curve from this observation and 
found marginal evidence for count rate variability (see 
Figure Set 13). Given the large number of counts 
available, we split the observation into 5 ks-long sub- 
exposures and extracted spectra from each sub-exposure 
to search for evidence of underlying spectroscopic vari- 
ability. Throughout the observation, a two-temperature 
thermal plasma model was sufficient to describe the X- 
ray spectrum of 6 Com, with the best fit temperatures 
and normalizations remaining constant within the un- 
certainties. All best-fit sub-exposure spectral parame- 
ters for 6 Com are provided in Table 5, and the average 
best-fit spectral parameters are listed in Table 6. The 
best-fit spectral models are shown in Figure Set 14. 


A.13. LHS 208 (x Men) 


All spectral parameters for LHS 208 are taken from 
Huang et al. (2018). LHS 208 hosts three confirmed 
exoplanets. HD 390931 b is a giant planet that was 
detected via RV measurements (Jones et al. 2002); it 
has Mp sini~12 Mjyp and an orbital period of ~2090 
days (corresponding to an orbital distance of ~3.3 AU; 
Feng et al. 2022). m Men c is a super-Earth that was 


discovered via transits (Gandolfi et al. 2018); it has 
an M,sini~3.5 Mẹ with an orbital period of ~6 days 
(~0.07 AU orbital distance; Feng et al. 2022). m Men 
d was detected via RV measurements with Mp sini~13 
Me and an orbital period of ~125 days (Hatzes et al. 
2022). LHS 208 not detected in a 36.9 ks observation 
with XMM-Newton or a 19.6 ks observation with Chan- 
dra. 


A.14. p CrB 


Some minor discrepancies, mass ranges from 0.95 
(Brewer et al. 2023) to 1.05 Mo (Stassun et al. 2017) 
and spectral type from GOV (Gray et al. 2006) to G2V 
(Turnbull 2015). p CrB hosts three confirmed exoplan- 
ets (p CrB c, d, and e), with the existence of a fourth 
exoplanet controversial (p CrB b, Brewer et al. 2023, 
which may simply be the result of correlated stellar ac- 
tivity in the RV measurements;). The masses of p CrB c, 
d, and e are Mp sin? ~28 Mg, ~22 Mg, and 3.8 Ms, re- 
spectively, with orbital periods of ~102 days (70.4 AU), 
7282 days (~0.8 AU), and ~13 days (~0.11 AU), re- 
spectively. p CrB is also approaching the end of its main 
sequence lifetime, with a prognosis of planetary engulf- 
ment for most of the known planets in the system (Kane 
2023). p CrB was observed once with Chandra/ ACIS-S 
for 9.8 ks, but was not detected. 


A.15. GL 672 


There are no significant discrepancies in the stellar 
physical parameters for GL 672, and there are no known 
or candidate exoplanets in the system. GL 672 was ob- 
served for 9.9 ks by Chandra/ ACIS-S but was not de- 
tected. 


A.16. x! Ori 


xl Ori is a GOV-type star in a ^14 year elliptical 
(e ~ 0.45) orbit with a low-mass companion. RV mea- 
surements compiled by Han & Gatewood (2002) were 
used to infer a secondary mass of ~0.15 Mo with a 
likely spectral type of MAV or M5V. The estimated an- 
gular semimajor axis of the orbit is ~0.668”, unresolv- 
able by XMM-Newton but potentially spatially resolv- 
able by Chandra (see Section 3). The inferred spectral 
parameters may be biased due to the presence of this 
low-mass companion. There is significant disagreement 
in the literature about the age of the system. Turnbull 
(2015) and Takeda et al. (2007) estimate an age of 4.32 
Gyr, while Mamajek & Hillenbrand (2008) use activity- 
rotation diagnostics to estimate an age of 300-400 Myr. 
'There are no confirmed or candidate exoplanets in the 
system. 

Since only one XMM-Neuton observation is available 
for x! Ori, we are unable to identify the source of the 
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X-ray emission — whether it be from the GOV star, the 
low-mass companion, or both. There is one potential, 
minor flaring event during the observation, starting ~21 
ks after the beginning of the observation and lasting for 
~4 ks, where the peak count rate exceeds the median 
count rate by ~40% (see Figure Set 13). While this may 
not be a bona fide flare event, we classify this period as 
flaring so that its spectrum can be compared to the re- 
mainder of the observation. The light curve otherwise 
shows a count rate that is either constant within the un- 
certainties (which we classify as quiescent) or somewhat 
elevated above the median count rate level by ~20% 
(which we classify as elevated). 

Figure Set 14 shows the quiescent, elevated, and po- 
tential flare spectra for the x! Ori system. The el- 
evated and flare spectra are nearly indistinguishable, 
and both show excess emission over the quiescent spec- 
trum at wavelengths <10A. This difference in short- 
wavelength emission is driven primarily by the hottest 
thermal plasma component in our model, which has a 
temperature of ~0.8 keV for the elevated and flaring 
spectra but ~0.7 keV in the quiescent spectrum (see 
also Table 6). The integrated X-ray luminosity of x! 
Ori is constant within the uncertainties across the three 
spectra. 


A.17. GL 788 


There are no significant discrepancies in reported stel- 
lar parameters of GL 788. The system does not host any 
known or candidate exoplanets. GL 788 was detected 
in a 9.6 ks XMM-Newton observation with ~4470 net 
counts. There was no evidence of count rate variability 
apparent in the X-ray light curve (see Figure Set 13). 
The spectrum was well-described with a single, cool 
thermal plasma (KT ~0.2 keV). 


A.18. 47 UMa 


There are no significant discrepancies in the stellar 
physical parameters reported in the literature for 47 
UMa. The star hosts three confirmed giant exoplanets 
that were detected via radial velocity measurements. 47 
UMa b, c, and d have masses of Mp sin i~2.4 Myup, ~0.5 
Myup, and 1.5 ~2.4 Myjup, respectively, and have or- 
bital periods (distances) of ~1077 days (~2 AU), ~2287 
days (77.8 AU), and —19000 days (712.4 AU), respec- 
tively. All planetary parameters were taken from Rosen- 
thal et al. (2021). 

47 UMa was detected at ~12o0 significance in a 6.2 ks 
XMM-Newton observation. We extracted a light curve 
from the observation but found no evidence of variability 
(see Figure Set 13), and the star was too faint to enable 
spectroscopic modeling. 
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A.19. GL 620.1A 


GL 620.14 is either co-moving with (Mugrauer 2019) 
or in a wide binary with (Holberg et al. 2002) the DA- 
type white dwarf WD 1620-391. The angular distance 
between the two components is ~5.8’, which would cor- 
respond to a physical separation of ~4500 AU at the 
distance of GL 620.1A. GL 620.14 is classified as a bar- 
ium star, and its over-abundance in s-process elements 
has been attributed to pollution due mass transfer dur- 
ing the white dwarf progenitor's AGB phase (Porto de 
Mello & da Silva 1997). There is some disagreement 
in the literature over the spectral type and age of GL 
620.1A. Spectral types range from G1V (Gray et al. 
2006) to a later G3/5V (Turnbull 2015); we adopt a 
spectral type G1V. While Turnbull (2015) quote an age 
of 2 Gyr, other estimates of the age of GL 620.1A are 
lower: Ghezzi et al. (2010) and Takeda et al. (2007) 
derive upper limits of «1 Gyr and «0.68 Gyr, respec- 
tively, while Sanz-Forcada et al. (2010) adopt an age of 
0.4 Gyr. One giant exoplanet is confirmed via RV mea- 
surements of GL 620.1A (Mayor et al. 2004); the planet 
has M,sini~1.2 Mjyp and an orbital period of ~528 
days (orbital distance of ~1.3 AU). 

GL 620.1A is detected in both publicly available 
XMM-Newton observations (Figure 12) with sufficient 
counts to enable spectral modeling (see Figure Set 14 
and Table 6). There is no significant variability present 
in the count rate light curves of GL 620.1A (Table 4, 
see Figure Set 13), so we do not attempt to divide the 
observations into sub-exposures. 


A.19.1. The White Dwarf WD 1620-391 


While GL 620.1A is robustly detected in both pub- 
licly available XMM-Newton observations, WD 1620- 
391 is not (see Figure 12a). The 30 upper lim- 
its on the count rates of 0.014 ct s^! (observation 
ID 0822070201) and 0.0236 ct s^! (observation ID 
0822070301). We use WebPIMMs’° to convert the ob- 
served XMM-Newton/PN count rates (both observa- 
tions were taken with the thick filter) to a flux assuming 
an APEC thermal plasma model with kT z 0.5 keV (the 
approximate best-fit model for the accreting white dwarf 
G 29-38, which is detected in X-rays by both XMM- 
Newton and Chandra; Cunningham et al. 2022; Estrada- 
Dorado et al. 2023). We find flux (luminosity) upper 
limits of 2.0x 10714 erg s~! cm-? (4x10?6 erg s71) and 
3.4x 10714 erg s71 cm? (6.8x 107 erg s71) for observa- 
tions 0822070201 and 0822070301, respectively. These 
luminosity limits are a factor of ~5-7 higher than the 


10 See https: //heasarc.gsfc.nasa.gov/cgi- bin/Tools/w3pimms/ 
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observed luminosity of G 29-38; targeted observations 
with a more sensitive X-ray instrument are needed to 
determine if WD 1620-391 is emitting X-rays at a levels 
comparable to other X-ray-detected white dwarfs. 


A.20. GL 311 (x! UMa) 


'There are no significant discrepancies in the stellar 
physical parameters reported in the literature for GL 
311. The system does not host any currently known or 
candidate exoplanets. 

The variability metrics for GL 311 (see Table 4) 
suggest that count rate variability is observed in the 
available XMM-Newton observation (observation num- 
ber 0111400101; see Figure Set 13). There are two mi- 
nor flare-like events in the light curve. The first begins 
at ~21 ks after the beginning of the observation, and 
reaches a peak count rate that is a factor of ~2 higher 
than the median count rate. 'This event lasts for ap- 
proximately 8 ks and it immediately followed by a sec- 
ond count rate increase, which exceeds the median count 
rate by a factor of ~1.6 and lasts for approximately 5 ks. 
The post-flare count rate is elevated by ~15% relative to 
the pre-flare count rate. Spectroscopic modeling before, 
during, and after the flaring event show that the tem- 
peratures of all three thermal components increased in 
temperature during the flaring event, although a modest 
decrease in the normalization (related to the VEM) dur- 
ing the flare led to a constant (within the uncertainties) 
luminosity (Figure Set 14; see also Table 6). 

'Thus, although GL 311 clearly exhibits a degree of 
rapid X-ray variability (on a timescale of a few hours), 
this activity does not have a significant affect on the 
overall X-ray luminosity output of the star. The relative 
stability of the star’s X-ray luminosity is also apparent 
in Figure 9, where the luminosity curve for GL 311 is 
near-vertical. 


A.21. a Cen 


a Centauri A (G2V) and B (K1V) are members of 
the extensively studied triple system that also includes 
Proxima Centauri (M5.5V; see Section A.56). a Cen A 
and B are in a ~80 year elliptical orbit (e~0.5; Akeson 
et al. 2021). One candidate exoplanet, œ Cen B b, was 
reported (Dumusque et al. 2012) but was later found to 
be a false positive signal (Rajpaul et al. 2016). There 
are currently no other known or candidate exoplanets in 
the a Cen AB system. 

The semi-major axis of the a Cen system is ~18”, but 
due to the brightness of the stars there is considerable 
blending in the XMM-Newton images (see Figure 12). 
We therefore do not attempt to extract X-ray products 
(light curves, spectra) for each star individually. The 


a Cen system exhibits significant X-ray variability (see 
Figure Set 13). The system is well-described by a three- 
component thermal plasma model, and increases in the 
observed count rates are driven largely by increases in 
the hottest plasma temperature (from ~0.6 keV during 
quiescent periods to ~0.8-0.9 keV during times associ- 
ated with flaring events) and normalization (which in- 
creases by a factor of ~12). The spectra of a Cen are 
shown in Figure Set 14 (see also Table 6). Due to the 
significant blending of the two components, we do not 
use the quiescent œ Cen spectrum as a baseline for in- 
ferring other stellar luminosities. 


A.22. 44 Boo 


44 Boo (also called ; Boo) is a triple star system with a 
GOV primary (44 Boo A) in ~210 year orbit with a spec- 
troscopic binary composed of a K0V (44 Boo Ba) dwarf 
and a K4V (44 Boo Bb) dwarf (Zasche et al. 2009; Zirm 
2011). 44 Boo A is separated from the spectroscopic 
binary by —0.2". The 44 Boo B binary has an orbital 
period of only —6.4 hours and is a possible contact bi- 
nary (Lu et al. 2001). There are no known or candidate 
exoplanets in the 44 Boo system. 

44 Boo is robustly detected with XMM-Newton in an 
18.8 ks observation, with count rate variations clearly 
seen in the light curve (see Figure Set 13). We do not 
include the final ~5 ks in our spectral analysis due to 
a significant background flaring event during the obser- 
vation. Our spectral modeling revealed no significant 
differences between the minor flaring event from ~6-10 
ks and the relatively constant count rate period from 
15-20 ks. The spectrum of 44 Boo did change during 
periods of declining count rate (the first 6 ks and from 
10-15 ks) as the result changes in the hottest thermal 
plasma component: k73 decreased from ~2.1-2.2 keV 
during quiescent/flaring times to —1.5 keV during pe- 
riods of secular decline, while the normalization of this 
component (related to the VEM) more than doubled. 
These changes to APEC #3 result in a decrease in X- 
ray flux at the shortest wavelengths (<10 A); see Table 6 
and Figure Set 14. 


A.23. 18 Sco 


There are minor differences in the reported spectral 
type of 18 Sco: Turnbull (2015) quote a spectral type of 
G5V, but Gray et al. (2006) classifies the star as a G2V. 
We adopt a spectral type of G2V. 18 Sco hosts one candi- 
date super-Earth with Mj sin? ~6.8 Mg and an orbital 
period of ~20 days (Laliotis et al. 2023). The star is de- 
tected in both a 20.1 ks XMM-Newton observation and 
a 62.2 ks observation with Chandra/ ACIS-S. However, 
the 18 Sco is not bright enough in either observation for 
timing or spectroscopic analysis to be performed. 
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A.24. 51 Peg 


There are no significant discrepancies reported for the 
stellar parameters of 51 Peg. The star is famous for 
hosting the first discovered exoplanet (Mayor & Queloz 
1995), the hot Jupiter 51 Peg b with Mp sin i ~0.46 Mjup 
and an orbital period of ~4.2 days (70.05 AU Rosen- 
thal et al. 2021). 51 Peg was observed with both Chan- 
dra and XMM-Newton, but was not detected in either 
observation. 


A.25. C Ret 


The wide binary ¢ Ret contains Ç! Ret, which is listed 
as a G3V/G5V by Turnbull (2015) but G2V as Gray 
et al. (2006), and ¢? Ret, a G2V star (del Peloso et al. 
2000). The two stars are separated by ~5.2’. The sys- 
tem does not hot any known or candidate exoplanets. 

Although this separation would be easily resolvable by 
XMM-Newton, only ¢! Ret is contained with the fields 
of view of the XMM-Newton cameras. The star is ro- 
bustly detected (with ~7700 net counts) in a 20.9 ks 
XMM-Newton observation (ID 0404920101). No signif- 
icant variability is observed in the X-ray light curve of 
Ç! Ret (see Figure Set 13), and the spectrum is well- 
described by a two-component thermal plasma model 
(with temperatures of ~0.1 keV and ~0.4 keV; see Fig- 
ure Set 14 and Table 6). 


A.26. HD 136352 


'There are no significant discrepancies in the stellar 
physical parameters reported in the literature for HD 
136352. The star hosts three confirmed exoplanets, all 
of which have masses in the super-Earth/sub-Neptune 
regime. HD 136352 b, c, and d have Mj, sini of ~4.7 
Ma, ^11 Mẹ, and ~8.8 Ma, respectively, with orbital 
periods (distances) of ~11.6 days (~0.1 AU), ~28 days 
(~0.17 AU), and ~107 days (70.43 AU), respectively. 
The inner two planets were found to transit the host star 
by Kane et al. (2020b), providing bulk densities for the 
planets. Subsequent observations by Delrez et al. (2021) 
revealed that the third, outer planet also transits the 
host star. HD 136352 is marginally detected (—2.76) in 
a 16.9 ks XMM-Newton observation with a count rate of 
0.0041+0.0015 ct s^! with the PN camera and medium 
filter (Table 9). 


A.27. GL 327 


'There are no significant discrepancies in reported stel- 
lar parameters of GL 327 in the literature. The system 
does not host any currently known or candidate exo- 
planets. GL 327 was robustly detected with ~10,500 
net counts in a 23.5 ks exposure with XMM-Newton 
(observaton ID 0404920201). Count rate variability 
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was observed in the X-ray spectrum as shown in Fig- 
ure Set 13, so we split the observation into three sub- 
exposures and modeled the three spectra independently. 
In all three sub-exposures, a single-temperature ther- 
mal plasma component was sufficient to describe the 
observations. During the quiescent period from ~6-21 
ks, the plasma temperature was ~0.4 keV, which in- 
creased to ~0.5 keV during the descending and flaring 
periods. The normalization during periods of increased 
X-ray count rate additionally increased by ~20%. The 
spectra of GL 327 are shown in Figure Set 14 (see also 
Table 6. 


A.28. u Ara 


No significant discrepancies in stellar parameters are 
present in the literature for u Ara, although there is 
some uncertainty about the evolutionary stage; we adopt 
a spectral type of G3IV/V (Turnbull 2015). u Ara hosts 
four known (giant) exoplanets, all detected via RV mea- 
surements. HD 160691 b, c, d, and e have masses of 
Mp sini~1.7 Mjup, ~2.4 Myup, ~10.2 Me, and ~0.5 
M5up, respectively, and the planets have orbital peri- 
ods (orbital distances) of ~646 days (~1.5 AU), ~4470 
days (~5.5 AU), ~10 days (~0.09 AU), and ~307 days 
(~0.9 AU), respectively. All planetary data were taken 
from GoZdziewski et al. (2007). 4 Ara is not detected 
in either of the two publicly available XMM-Newton ob- 
servations. 


A.29. k! Cet 


'There is some discrepancy in the reported age esti- 
mates of x! Cet: Turnbull (2015) lists an age of 2.2 Gyr 
and Takeda et al. (2007) derive an upper limit of «2.76 
Gyr, but many studies report significantly younger ages. 
For example, Mamajek & Hillenbrand (2008) derive an 
age of 300-400 Myr, while Güdel et al. (1997) estimated 
an age of 750 Myr and and Dorren & Guinan (1994) es- 
timate an age of 650 Myr. There are no currently known 
or candidate exoplanets orbiting &! Cet. 

«! Cet is relatively well-observed in X-rays, with five 
XMM-Newton observations and two Chandra observa- 
tions available in the archives, and the star is easily de- 
tected in all exposures. It was also observed twice with 
Chandra using ACIS-S7 in 1/8 subarray mode with a 
0.44 s readout time in both observations. While there 
are clear indications of variability in the light curves of 
k! Cet the maximum observed count rate only ever ex- 
ceeds the median count rate by a factor of —2 (see Fig- 
ure Set 13). However, despite the only modest changes 
in count rate, the spectra during these minor flare-like 
events is quite distinct from the lower count rate qui- 
escent periods, with significant excess emission a wave- 
lengths shortward of 10 À compared to quiescent times. 


The change in count rate is being driven by a signifi- 
cant increase in the temperature (by ~0.2 keV, or ~2.3 
MK) and normalization (a factor of —2.2 increase) of 
the hottest component in a three-thermal plasma com- 
ponent model (Figure Set 14; see also Table 6). 


A.30. HD 140901 


There are no significant discrepancies in the published 
stellar physical parameters of HD 140901. Two planets 
have been identified via RV measurements. HD 140901 c 
has Mp sin i ~1.8 Mjup with an orbital period of ~14390 
days (orbital distance of ~11.8 AU; Philipot et al. 2023). 
The existence of HD 140901 b is controversial; if the 
planet exists, it has a mass of My sini of ~16Mẹ with 
an orbital period of ~9 days (orbital distance 0.085 AU), 
but the RV data of HD 140901 can be equally well 
described by a retrograde orbital inclination of ~138° 
(Feng et al. 2022; Philipot et al. 2023). 

HD 140901 was observed once by Chandra/ ACIS-I for 
24.7 ks. The star was robustly detected with ~220 net 
counts. The light curve shows clear evidence for vari- 
ability (see Figure Set 13), but the star is not bright 
enough to enable spectral modeling. 


A.31. GJ 777 


GJ 777 is a binary star system composed of a G6IV 
primary (GJ 777A) and an M4.5V secondary (LHS 
2509), which are separated by ~3000 AU (~3’). GJ 
777A hosts two currently known exoplanets, detected 
via RV measurements: HD 190360 b has My sini ~1.8 
Myup and an orbital period of —2854 days (~3.9 AU; 
Feng et al. 2021), and HD 190360 c has M, sin i ~19 Mg 
and an orbital period of ~17 days (70.13 AU; Rosenthal 
et al. 2021). No X-ray emission is detected coincident 
with the GJ 777 binary system in either of two snapshot 
XMM-Newton observations. 


A.32. € Boo 


€ Boo is a visual binary (€ Boo A has a spectral type 
G8V and € Boo B has a spectral type K4V; Turnbull 
2015) with an angular separation of ~4.9” (orbital pe- 
riod ~150 years; Wielen 1962), the stars are unresolvable 
by XMM-Newton (see Figure 12). No discrepancies in 
stellar physical parameters in the literature. No known 
or candidate exoplanets. € Boo A is a known BY Dra- 
conis variable (Finley et al. 2019). 

The € Boo binary is robustly detected by XMM- 
Newton, although the two components cannot be re- 
solved. The system exhibits dramatic X-ray variability, 
as shown in the light curve in Figure Set 13. The X- 
ray spectra of € Boo are described by three-component 
thermal plasma models in all variability states. The flar- 
ing spectrum requires hotter temperatures for all three 


plasma components than the quiescent spectrum, with 
the normalization of the hottest component increasing 
by a factor of ~3.9 during flaring times compared to qui- 
escent times. Figure Set 14 shows the best fit spectral 
models for € Boo (see also Table 6). 


A.33. 6 Pav 


There are no significant discrepancies in the stellar 
physical parameters reported in the literature for 6 Pav. 
The star does not host any currently known or candidate 
exoplanets. 

ó Pav was detected in a 32.6 ks XMM-Newton obser- 
vation (ID number 0780510401) ~1000 net counts. The 
A? statistic provides evidence that the star may be vari- 
able on short timescales, however there are not a suffi- 
cient number of counts to split the observation into sub- 
exposures (see Figure Set 13). We attempted to model 
the spectrum of 6 Pav despite the low number of counts; 
the spectrum required two thermal plasma components 
to achieve a statistically acceptable fit (see Table 6). 
The best-fit spectrum in shown in Figure Set 14. 


A.34. LHS 2156 


There are no significant discrepancies in the stellar 
physical parameters reported in the literature for LHS 
2156. There are no known or candidate exoplanets in 
the system. 

LHS 2156 was detected with ~60,000 net counts in 
a 68.4 ks XMM-Newton observation with the PN cam- 
era (thick filter). Figure Set 13 shows the X-ray light 
curve, which exhibits count rate variability over the du- 
ration of the observation. The spectra of LHS 2156 are 
well described by three-component thermal plasma mod- 
els. The best-fit temperatures and normalizations of 
two cooler plasma components remain constant within 
the uncertainties in all sub-exposures. The hottest ther- 
mal plasma component increases in temperature by ~0.1 
keV during flaring periods compared to quiescent peri- 
ods, an the normalization increases by ~25%. The best- 
fit spectra for LHS 2156 are shown in Figure Set 14 (see 
also Table 6). 


A.35. 82 Eri 


There is some discrepancy in the literature about 
the age of 82 Eri: Turnbull (2015) and Takeda et al. 
(2007) quote ages >12 Gyr, but a much younger age of 
5.7640.66 Gyr was found by (Pepe et al. 2011). 82 Eri 
hosts four confirmed exoplanets, all of which were dis- 
covered via RV measurements (Pepe et al. 2011). HD 
20794 b, c, d, and e have My sin? values of ~2.8 Mẹ, 
72.5 Mg, ~3.5 Ma, and ~4.8 Ma, respectively, with or- 
bital periods (distances) of ~18 days (~0.13 AU), ~43 
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days (~0.23 AU), ~89 days (~0.26 AU), and ~147 days 
(~0.5 AU), respectively. All planetary parameters were 
taken from Feng et al. (2017a). 82 Eri has been observed 
by both XMM-Newton and Chandra (see Table 3), but 
was not detected in either observation. 


A.36. 7 Ceti 


There are no significant discrepancies in the stellar 
physical parameters of 7 Ceti reported in the litera- 
ture. 7 Ceti hosts four exoplanets, all confirmed via 
radial velocity measurements (Feng et al. 2017b), with 
some evidence of additional planets in or near the r Cet 
HZ (Dietrich & Apai 2021). 7 Ceti e and f both have 
masses Mysini ~3.9 Ma with orbital periods of ~163 
days (~0.55 AU) and ~636 days (~1.3 AU), respec- 
tively, while r Ceti g and h both have masses of Mp sin i 
71.8 Mg and orbital periods of ~20 days (~0.13 AU) 
and ~49 days (~0.24 AU), respectively. 

T Ceti was observed by Chandra ACIS-S for 45.1 ks 
in 1/8 subarray mode and detected at ~7o significance 
with ~110 net counts. We extracted a light curve from 
this observation and find an A? statistic that is strongly 
suggestive of variability (see Table 4 and Figure Set 13), 
but due to the relatively low number of counts we did not 
attempt to extract a spectrum. 7 Ceti was marginally 
detected in an 11.9 ks XMM-Newton observation, but 
with insufficient counts for light curve extraction or 
spectral modeling to be performed. The luminosity esti- 
mates for 7 Ceti suggest that the star’s X-ray luminosity 
changes by a factor of ~2. 


A.37. 55 Cnc 


55 Cnc is a binary star system containing a G8V star 
(55 Cnc A) and a M4.5V dwarf (55 Cnc B). No signif- 
icant discrepancies in stellar physical parameters. The 
two components are separated by ~1,065 AU (~1.4’; 
Raghavan et al. 2006). Five exoplanets are known to 
orbit 55 Cnc A. 55 Cnc b was the first to be discovered 
(Butler et al. 1997) with a mass of My sini ~0.8 Myup 
and an orbital period ~14.7 days (corresponding to an 
orbital distance of ~0.11 AU; Bourrier et al. 2018; Nel- 
son et al. 2014). 55 Cnc c, d, e, and f has masses of 
Mypsini ~0.16 Mjup, 73.9 Mjup, ~8 Me, and ~0.15 
M;jyp, respectively, with orbital periods of ~44 days 
(~0.24 AU), ~4870 days (~5.5 AU), ~0.7 days (~0.02 
AU), and ~260 days (~0.8 AU), respectively. 

Cnc 55 was observed twice by Chandra/ACIS-S in 1/8- 
subarry mode (for 10.7 ks and 18 ks) and once for 10.4 ks 
by XMM-Newton. 55 Cnc A was not detected by Chan- 
dra, but was detected (at ~5.6c) significance by XMM- 
Newton. No evidence for strong variability was observed 
in the XMM-Newton light curve of 55 Cnc A (Table 4, 
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see also Figure Set 13). Despite the relatively low num- 
ber of net counts, we attempted to model the X-ray 
spectrum of 55 Cnc A. We find that a single-component 
thermal plasma model is sufficient to describe the X-ray 
spectrum (see Table 6 and Figure Set 14). 55 Cnc B was 
not detected in the XMM-Newton observation or one of 
the Chandra observations (the star was not within the 
field of view of the second Chandra observation). 


A.38. 70 Oph 


70 Oph is a visual and spectroscopic binary containing 
a KOV (70 Oph A) star and a K5V (70 Oph B) star with 
an orbital period of ~88 years (Pourbaix 2000). The or- 
bit is highly elliptical, with an eccentricity e = 0.50 that 
causes the orbital separation of the stars to vary from 
~11.4 AU to ~34.8 AU (Eggenberger et al. 2008). These 
orbital separations correspond to angular separations 
of ~2.3-6.8”, unresolvable by XMM-Newton but poten- 
tially spatially resolvable by Chandra (see Section 3). 
Eggenberger et al. (2008) find a best-fit age of 6.2+1.0 
Gyr, in broad agreement with the 6.8-7.0 Gyr found by 
Tang et al. (2008). These ages (which where derived 
using astroseismic observations) are in contrast to the 
significantly younger age (~680 Myr) predicted by Ma- 
majek & Hillenbrand (2008). There are no currently 
known or candidate exoplanets in the system. 

There are three XMM-Newton observations available 
for the 70 Oph system. The individual stars are not 
resolvable in the XMM-Newton observations (see Fig- 
ure 12); we therefore cannot separate the X-ray emission 
from the KOV star (70 Oph A) from the K5V star (70 
Oph B) with the currently available observations. 

Despite the changes in X-ray count rate observed in 
the light curves (see Figure Set 13), the best-fit spectral 
models of 70 Oph across different VarFlags are nearly 
indistinguishable (Figure Set 14; see also Table 6), and 
the X-ray luminosity of the two stars changes only by 
~40% between the flaring and quiescent periods. Since 
we cannot separate the X-ray emission of the binary 
components individually, we do not use the best-fit spec- 
trum of the 70 Oph system as template for inferring 
luminosities of other stars in our sample. 


A.39. 40 Eri 


The 40 Eri triple system contains a KOV star (40 Eri 
A) that is either in a wide orbit or co-moving with a 
white dwarf (40 Eri B)-M4.5V dwarf (40 Eri C) pair. 
The projected angular distance between 40 Eri A and 40 
Eri BC is ~1.3’. While the orbital period of the wider 
A-BC pair is very long, the orbital period of the B-C pair 
is measured to be ~230 years (Mason et al. 2017; Bond 
et al. 2017). The angular separation between 40 Eri B 


and 40 Eri C is ~6.9”, easily resolvable by Chandra (see 
Figure 2). 40 Eri C has a visual magnitude of ~11, which 
we use to estimate Ly, and Johnson & Wright (1983) 
estimate Tig to be ~3,100 K. There is one confirmed 
planet orbiting 40 Eri A, which was discovered via RV 
measurements. HD 26965 b has M,sinz of ~8.5 Ma 
with an orbital period of ~42 days (Ma et al. 2018). 
No currently known or candidate planets are associated 
with the 40 Eri BC pair. 

The first soft (0.1-2.0 keV) X-ray detection of the 40 
Eri system was obtained by HEAO-1, which measured a 
luminosity of (9.6--3.2) x 1078 erg s^! (corrected for dif- 
ferent distance assumed in that work; Cash et al. 1979). 
Due to the poor spatial resolution of HEAO-1, the bi- 
nary components could not be resolved. Different sce- 
narios for X-ray production are explored in Cash et al. 
(1979), and the X-ray emission was attributed to 40 Eri 
C. 

'The Chandra image of the system clearly shows X-ray 
sources coincident with both 40 Eri A and C (although 
not the white dwarf, 40 Eri B; see discussion below). 40 
Eri C is the brighter of the two sources, with roughly 
an order of magnitude more net counts than observed 
for 40 Eri A. There is no evidence for strong variability 
in either star (see Figure Set 13). While the coolest 
thermal plasma component in the 40 Eri C spectrum 
is similar in temperature to other late-type stars in our 
sample, the hotter two thermal plasmas are considerably 
hotter (see Table 6 and Figure Set 14). Despite the 
lack of obvious variability in 40 Eri C, the best-fit X-ray 
spectral parameters most closely resemble those of WX 
UMa in its flaring and post-flaring states. It is possible 
this star indeed exhibits significant coronal variability, 
but it was not detectable in the short snapshot Chandra 
observation available in the archive. 


A.39.1. The White Dwarf 40 Eri B 


'The white dwarf 40 Eri B is not detected in the 5 ks 
Chandra/ ACIS-S observation of the system. The 3e up- 
per limit on the count rate is 0.0033 ct s~'. We again 
use WebPIMMs to convert the observed count rate up- 
per limit to a flux (and luminosity) upper limit, assum- 
ing the same spectral model as for WD 1620-391 (Sec- 
tion A.19.1). We find a flux upper limit of 1.6x10 7? 
erg s^! cm"?, which corresponds to a luminosity of 
~5x10?5 erg s-!. This luminosity upper limit is be- 
low the measured X-ray luminosity of G 29-38 (78x 1075 
erg s !; Cunningham et al. 2022; Estrada-Dorado et al. 
2023), strongly suggesting a lack of a cold accretion disk 
around this white dwarf. 


A.40. 6 Eri 


There are no significant discrepancies in the reported 
physical parameters for 6 Eri, and the star does not host 
any known or candidate exoplanets. 6 Eri was observed 
once with XMM-Newton (observation observation) for 
56.2 ks. There was marginal evidence for significant 
variability in the X-ray light curve of ô Eri (see Figure 
Set 13), but the star was too faint for spectral model- 
ing. 


A.41. GL 451A (Groombridge 1830) 


There are no significant discrepancies in the stellar 
physical parameters reported in the literature for GL 
451A. The system does not host any currently known 
or candidate exoplanets. GL 451A was observed once 
with Chandra/ ACIS-S (observation ID 9931) for 32.8 ks 
in 1/8-subarray mode across ACIS detectors 5-7 with a 
readout time of 0.54 s. No X-ray source is found coinci- 
dent with the star's position. 


A.42. GL 117 


Turnbull (2015) report an age of 1.5 Gyr for GL 117 
and Takeda et al. (2007) derive an upper limit of «1.2 
Gyr. These ages are in contrast with the very young 
age (~100 Myr) predicted by Mamajek & Hillenbrand 
(2008). There are no currently known or candidate ex- 
oplanets orbiting GL 117. 

GL 117 is robustly detected in XMM-Newton observa- 
tion 0203060501, and there are low-level count rate vari- 
ations observable in the light curve (see Figure Set 13). 
The largest of these variations reaches a count rate level 
that is ~30% above the median count rate. Given 
the ample number of X-ray counts recorded by XMM- 
Newton, we divide the light curve into sub-exposures but 
we find no significant differences in the best fit spectral 
parameters nor the predicted luminosities (see Figure 
Set 14 and Table 6). The bright X-ray luminosity of 
the star, particularly when compared to other K dwarfs 
in our sample, are consistent with a young age for GL 
117. 


A.43. GL 788 


GL 783 is a binary system containing a K3V star (A) 
and an M4.5V star (B; Turnbull & Tarter 2003). The 
two stars were stated to be separated by ~8” by Allen 
et al. (2000) but Gaia DR3 positions indicate a sepa- 
ration of ~2.6” (Gaia Collaboration et al. 2023), both 
of which are unresolvable by XMM-Newton but poten- 
tially spatially resolvable by Chandra (see Section 3). 
GL 783B has an apparent visual magnitude of ~12.5 
(Allen et al. 2000), which we use to estimate Lyg. At 
the time of this writing there were no published values 
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of Teg for GL 783B. The system does not host any cur- 
rently known or candidate exoplanets. 

There is one XMM-Newton observation available for 
the GL 783 system, and the individual stars are not 
resolvable (see Figure 12). The therefore cannot sepa- 
rate the X-ray emission from the K3V star (GL 783A) 
from the M4.5V star (GL 783B) with the currently 
available observations. The A? statistic provides evi- 
dence that one or both stars may be variable on short 
timescales (see Figure Set 13), however there are not a 
sufficient number of counts to split the observation into 
sub-exposures. We attempted to model the spectrum of 
the GL 783 system despite the low number of counts; the 
spectrum required three thermal plasma components to 
achieve a statistically acceptable fit (see Figure Set 14 
and Table 6). The coolest thermal plasma component 
(~0.1 keV) is similar to the coolest components of the 
other FGK stars in our sample, while the temperatures 
of the hottest (^1 keV) and intermediate (70.3 keV) 
components more closely resemble components found in 
the later-K and M dwarfs in our sample. 

Since we cannot separate the X-ray emission of the bi- 
nary components individually, we do not use the best-fit 
spectrum of the GL 783 system as template for inferring 
luminosities of other stars in our sample. Follow-up ob- 
servations with higher spatial resolution are needed to 
associate the observed X-ray emission with one or both 
binary components. 


A.44. GL 183 


There is some discrepancy in the reported spectral 
type of GL 183: Gray et al. (2006) report a spectral 
type K3V, while K4III is reported by Turnbull (2015). 
We adopt a spectral type K3V. There is additional ten- 
sion in the age estimates, with a younger age (^2 Gyr) 
reported by Mamajek & Hillenbrand (2008) and older 
ages reported by Turnbull (2015, 5.3 Gyr) and Takeda 
et al. (2007, «5.45 Gyr). There are no currently known 
or candidate exoplanets orbiting GL 183. 

GL 183 was detected in a 18.4 ks XMM-Newton obser- 
vation (ID number 0780510301) ~650 net counts. The 
A? statistic provides evidence that the star may be vari- 
able on short timescales (see Figure Set 13), however 
there are not a sufficient number of counts to split the 
observation into sub-exposures. We attempted to model 
the spectrum of GL 183 despite the low number of counts 
and found the spectrum required two thermal plasma 
components to achieve a statistically acceptable fit (Fig- 
ure Set 14; see also Table 6), and the fit parameters were 
reasonably well-constrained. 
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A.45. GJ 667 


The GJ 667 is a triple star system wherein the two 
more massive components, A and B (both of which are 
early/mid K dwarfs) orbit one another in a highly ellip- 
tical orbit (e ~ 0.58). The average angular separation of 
the two stars is ~1.8”, possibly resolvable by Chandra. 
The third star in the system, the M1.5V dwarf GJ 667 C, 
is separated from the AB pair by ~0.5’ (which would be 
easily resolvable with both Chandra and XMM-Newton). 
There is no published age estimate for the system. GJ 
667C has been studied previously in X-rays by Brown 
et al. (2023), and is known to host two confirmed super- 
Earth exoplanets: GJ 667C b, with M, sini ~5.6 Ma 
and an orbital period of ~7.2 days (70.05 AU; Robert- 
son & Mahadevan 2014; Bonfils et al. 2013), and GJ 
667C c, with with Mj sini ~3.8 Ma and an orbital pe- 
riod of ~28 days (~0.125 AU; Anglada-Escudé et al. 
2013). An additional three planets in wider orbits were 
reported by Anglada-Escudé et al. (2013) but were un- 
able to be confirmed by Robertson & Mahadevan (2014); 
their existence remains controversial. 

GJ 667 was observed twice with Chandra/ ACIS-S in 
1/8-subarray mode. GJ 667 A and B were not within 
the field of view of either observation. GJ 667 C was 
detected in both observations. We extracted an X-ray 
light curve from Chandra ObsID 17317 (with an expo- 
sure time of 18.2 ks, see Figure Set 13), but the star was 
too faint for spectral modeling to be performed. There 
were insufficient counts in the second Chandra observa- 
tion (ObsID 17318, with an exposure time of 9.1 ks) for 
light curve or spectral analysis. 


A.46. LHS 1875 


'There is some discrepancy in the literature about the 
age of LHS 1875: Turnbull (2015) quote an age of ~3.1 
Gyr, while Takeda et al. (2007) provide a significantly 
younger upper limit of «0.48 Gyr. There are no known 
or candidate exoplanets orbiting LHS 1875. LHS 1875 
was detected with ~9260 net counts in a 21.4 ks XMM- 
Newton/PN observation. No evidence for significant X- 
ray variability was observed (see Figure Set 13), and 
the spectrum was well-described by a three-component 
thermal plasma model (see Table 6). 


A.47. GL 570 


There is a significant age discrepancy in the literature 
for the GL 570 quadruple-star system: Takeda et al. 
(2007) report a very young age of <0.6 Gyr, while Turn- 
bull (2015) adopt an ade of ~3 Gyr. There are no other 
significant discrepancies in the stellar physical parame- 
ters for the stars in this system. There are no currently 
known or candidate exoplanets in the system. 


GL 570A (K4V) orbits a binary red dwarf pair (GL 
570B, an M1.5V dwarf, and GL 570C, an M3V dwarf) 
with a semi-major axis of ~190 AU (corresponding to 
an angular separation of ~0.5’). GL 570A and the GL 
570BC pair are both detected in X-rays and are easily 
resolvable with Chandra (see Figure 12). The red dwarfs 
the make up the BC pair are separated by ~0.2” (Mar- 
iotti et al. 1990; Forveille et al. 1999) and cannot be 
resolved. The fourth component, GL 570D, is a brown 
dwarf in a very wide orbit (~4’) from the GL 570 ABC 
triple system (Burgasser et al. 2000). GL 570D is not 
detected in X-rays. 

We extract X-ray light curves for both GL 570A and 
the GL 570BC pair. There is clear evidence for variabil- 
ity in both light curves (see Figure Set 13), although 
there are not enough net counts from GL 570BC to 
enable spectroscopic modeling on sub-exposures. We 
therefore divide observation of GL 570A into three dis- 
tinct sub-exposures (the first ^20 ks, 20-30 ks, and 30-40 
ks) and independently model the spectra, but only ex- 
tract one spectrum of GL 570BC from the entire ob- 
servation. All spectra are well described by a two- 
component thermal plasma model. The flaring event 
seen from 20-30 ks in the GL 570A light curve is charac- 
terized by changes in plasma temperatures (70.1 keV for 
the cooler component and ~0.2 keV for the hotter com- 
ponent) and an increase in the hot-component normal- 
ization by ~170%. The different best-fit spectral models 
for GL 570A, as well as a comparison of the flaring-to- 
quiescent spectra, are shown in Figure Set 14. The GL 
570A model is overall hotter than the GL 570BC model 
(see Table 6). 


A.48. 61 Cyg 


'The 61 Cygni AB system is a wide binary containing 
a K5V star (A) and a K7V star (B), with an angular 
separation of ~12” (the orbital period of the system is 
7-659 years). Most age estimates suggest the system is 
old (641 Gyr; e.g., Kervella et al. 2008). The system is 
resolved by XMM-Newton (Robrade et al. 2012). The 
system does not host any currently known or candidate 
exoplanets. 

61 Cyg has been observed 38 times with XMM- 
Newton, making it the best X-ray-observed star system 
in our sample. The observations are spaced ~6 months 
apart and span ^18 years, and the two stars are well- 
resolved in all available images (see Figure 12). Both 
stars are observed to flare (see Figure Set 13), and all 
the spectra extracted from both stars are well describe 
by a three-component thermal plasma model. The three 
best-fit temperatures in the models of 61 Cyg A remain 
constant within the uncertainties during different vari- 


ability epochs. The changes in count rate observed dur- 
ing flaring events in 61 Cyg A are driven by changes 
in the normalizations of the two hotter thermal plasma 
components (and hence the VEMs): the hottest compo- 
nent normalization increases by an order of magnitude, 
and the intermediate temperature component normal- 
ization increases by a factor of ~2.5. The spectra of 
61 Cyg B show similar changes in normalization during 
flaring events, but the two cooler plasma components 
also increase in temperature during flaring events (by 
~0.1 keV). The best-fit spectra for 61 Cyg A and B are 
shown in Figure Set 14 (see also Table 6). 


A.49. GL 412 


The GL 412 AB system contains a M1V primary (A) 
and a UV Ceti-type flaring M6.6V secondary (B; Mann 
et al. 2015, also referred to as WX UMa). X-ray flares 
were previously observed from the system with ROSAT 
(Schmitt et al. 1995) and attributed to GL 412B; GL 
412A was not believed to be a significant source of X- 
rays. The system does not host any currently known or 
candidate exoplanets. 

The two stars have an angular separation of ~31” and 
are resolved in the XMM-Newton image (observation 
number 0742230101; see Figure 2). The image clearly 
suggests that the observation caught a significant WX 
UMa flaring event and that the M1V primary, while 
fainter, emits a detectable quantity of X-ray radiation. 
A third X-ray source is detected ~37” from WX UMa, 
but is not obviously associated with the GL 412 binary. 
We measure a count rate of 0.057+0.002 ct s^! in the 
PN image (medium filter) for GL 412A. 

The light curve of GL 412B is shown in Figure Set 13. 
Spectral modeling shows that the X-ray luminosity in- 
creases by a factor of ~30 between quiescent and flar- 
ing times. The quiescent spectrum (the first ~2 ks of 
the observation) is well-described with only two thermal 
plasma components. The steep, ~0.5 ks flare, the ele- 
vated count rate (from ~3-7 ks), and subsequent declin- 
ing period all require three components to adequately 
describe the X-ray spectra. The flaring spectrum ex- 
hibits a very hot thermal plasma component, with a 
temperature of ~3.5 keV (~40 MK), which cools as the 
flare subsides (it has a temperature of ~2.2 keV dur- 
ing the “elevated” period and ~1.4 keV during the “de- 
scending” count rate period; see Figure Set 14 and Ta- 
ble 6). 


A.50. GJ 832 


No significant discrepancies in stellar physical param- 
eters. One confirmed exoplanet, GJ 832 b with Mj sin 
70.6 Mjyp and an orbital period of 3853 days (corre- 
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sponding to an orbital semi-major axis of ~3.7 AU, Bai- 
ley et al. 2009; Philipot et al. 2023). 

GJ 832 was observed twice by XMM-Newton, once 
during a 8.9 ks (observation 0748010201) exposure dur- 
ing which the star was detected but too faint to en- 
able timing or spectroscopic modeling, and once during 
a 28.6 ks exposure (observation 0860303301) when the 
star was overall brighter and underwent a short flar- 
ing event (see Figure Set 13). The quiescent spectrum 
is well-described by two thermal plasma components, 
while the lower signal-to-noise during the short flaring 
period allows the flare spectrum to be modeled with a 
single thermal plasma component. Overall, the luminos- 
ity of GJ 832 increases by a factor of ~2.5 during the 
flare. The best-fit spectral models and a comparison of 
the flaring-to-quiescent spectrum are shown in Figure 
Set 14 (see also Table 6). The spectrum extracted dur- 
ing the short flaring time during the observation shows 
an excess of X-ray flux at ~2-30 A, but shows weaker 
X-ray emission (compared to the quiescent spectrum) at 
wavelengths <2 A. 


A.51. Kapteyn’s Star 


There are no significant discrepancies in the stel- 
lar physical parameters reported in the literature for 
Kapteyn’s Star. Although the detection of sub-Neptune 
mass planet was reported by Anglada-Escude et al. 
(2014), this was later refuted by Bortle et al. (2021) who 
argued that the observed RV variations were artifacts of 
stellar activity and rotation. There are not currently 
any additional candidate exoplanets around Kapteyn's 
Star. 

Kapteyn's Star has been observed numerous times by 
XMM-Newton, Chandra, and Swift (see Table 3), but is 
not detected in any of archival X-ray observations. We 
explored whether Kapteyn's Star would be detected in a 
stacked Chandra X-ray image (the Chandra background 
is lower and significantly less affected by spurious flaring 
events than XMM-Newton). The stacked Chandra image 
has an effective exposure time of 72.5 ks. Kapteyn's Star 
is still not detected in the stacked image, with a count 
rate upper limit of «1.24x10^^ ct s71. Assuming the 
best fit quiescent spectrum of Wolf 1055, this count rate 
corresponds to a luminosity upper limit of <1.1x10?4 
erg st. 


A.52. Wolf 1055 


There are no significant discrepancies in the stellar 
physical parameters of Wolf 1055 or its binary compan- 
ion VB 10 (GL 572B, a M8V-type BY Draconis flare 
star Burt et al. 2021). There is one confirmed exoplanet 
orbiting Wolf 1055: HD 180617 b (also identified as GL 
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572 Ab), with M, sini ~12 Mg and an orbital period of 
~106 days (~0.34 AU). 

Wolf 1055 and VB 10 are separated by ~1.2’ and are 
easily resolvable with both XMM-Newton and Chandra. 
Both Wolf 1055 was observed once with XMM-Newton 
(for 24.2 ks) and twice with Chandra (with exposure 
times of 12.2 ks and 29.2 ks). Both Wolf 1055 and VB 10 
are detected in all three observations. A comparison of 
the XMM-Newton observation and the longer Chandra 
observation (ObsID 7428) is shown in Figure 12. Wolf 
1055 is bright enough in two out of these three obser- 
vations to enable spectral modeling. Although VB 10 is 
bright enough to be significantly detected in all three ob- 
servations and we are able to extract light curves from 
all three observations. The light curves of both stars 
obtained during XMM-Newton observation 0504010101 
are shown in Figure Set 13. 

The quiescent spectrum of Wolf 1055 is well-described 
by a three-component thermal plasma model, while the 
lower signal-to-noise during the flaring period in the sec- 
ond half of the XMM-Newton observation 0504010101 
allows the flare spectrum to be modeled with a single 
thermal plasma component. Overall, the luminosity of 
Wolf 1055 increases by ~80% during the flare. The best- 
fit spectral models and a comparison of the flaring-to- 
quiescent spectrum are shown in Figure Set 14 (see also 
Table 6). The spectrum extracted during the flaring 
period during the second half of observation shows an 
excess of X-ray flux at ~2 A, and shows weaker X-ray 
emission (compared to the quiescent spectrum) at wave- 
lengths <2 A. 


A.53. AD Leo 


There are no significant discrepancies in the stellar 
physical parameters reported for AD Leo; all stellar pa- 
rameters are taken from Kossakowski et al. (2022) and 
references therein. There are no confirmed or candidate 
exoplanets orbiting AD Leo. 

AD Leo is a known flaring star, and exhibits strong 
X-ray flaring both times the star was observed by XMM- 
Newton. Both light curves are shown in Figure Set 13. 
High count rates enable spectroscopic modeling over 
multiple sub-exposures. The X-ray spectra during qui- 
escent periods are generally well-described with a two- 
component thermal plasma model, while periods of 
flaring and enhanced variability require three thermal 
plasma components to adequately describe the observed 
spectra. Increases in X-ray count rate are driven almost 
entirely by the appearance of this third, hot (~1.5 keV) 
thermal plasma component, which drives a significant 
increase (by three- to five-orders of magnitude) in X-ray 
emission below ~2 A (see Figure Set 14). 


A.54. Wolf 1061 


There are no significant discrepancies in stellar physi- 
cal parameters of Wolf 1061. No age estimate is available 
for the star. Wolf 1061 hosts three currently known exo- 
planets (Wright et al. 2016): Wolf 1061 b, c, and d have 
masses of Mp sin i ~1.9 Mg, ~3.4 Ma, and ~7.7 Mg, re- 
spectively, with orbital periods of ~4.9 days (70.04 AU), 
717.9 days (~0.09 AU), and ~217 days (~0.47 AU), 
respectively. All planetary data taken from Astudillo- 
Defru et al. (2017). 

Wolf 1061 was detected in one 38.3 ks observation with 
Chandra/ACIS-S with ~240 net counts. The light curve 
shows evidence for variability (see Figure Set 13), but 
the star is not bright enough to enable spectroscopic 
modeling. 


A.55. Luyten’s Star 


There are no significant discrepancies in the reported 
stellar parameters for Luyten’s Star in the literature. No 
age estimate is available. There are two confirmed exo- 
planets orbiting Luyten’s Star: GJ 273b and GJ 273c, 
which were detected via RV measurements by Astudillo- 
Defru et al. (2017). GJ 273b and GJ 273c have masses 
of ~2.9 Ma and ~1.2 Ms, respectively, and the planets 
have orbial periods of ~18.6 days (70.09 AU) and ~4.7 
days (70.03 AU), respectively. 

Luyten's star is detected at —5c significance in a 28.4 
ks observation with Chandra/ ACIS-S. The light curve 
shows evidence of X-ray variability (see Figure Set 13), 
but are are unable to extract X-ray spectra due to the 
low number of net counts. 


A.56. Proxima Centauri 


Proxima Centauri is a member of an extensively stud- 
ied triple system along with a Cen A and a Cen B. 
There is currently one confirmed Earth-like exoplanet, 
Proxima Centauri b, which has Mj sin? ~1 Ma and an 
orbital period of ~11 days (orbital distance ~0.05 AU; 
Faria et al. 2022; Damasso et al. 2020; Anglada-Escudé 
et al. 2016). Two candidate exoplanets, Proxima Cen- 
tauri c and d, are awaiting follow-up observations. If 
present, Prox Cen c (Damasso et al. 2020) and d (Faria 
et al. 2022) are estimated to have M, sini values of 
^6 Mẹ and ~0.26 Ma, respectively, with orbital pe- 
riods (distances) of ~1900 days (71.5 AU) and ~5 days 
(~0.029 AU). 

Proxima Centauri exhibits significant flaring events, 
such as those illustrated in Figure 3 (see also Figure 
Set 13). The X-ray spectrum of Prox Cen is well- 
described by a two-component thermal plasma model 
(with temperatures of — 0.25 keV and ~0.85 keV) during 
quiescent times and periods of elevated or rising count 
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Figure 15. Same as Figure 10 with Swift/XRT sample in 
light blue. 


rates. A third thermal plasma component is required to 
describe the spectra during flaring periods and during 
periods of descending count rates following a flare. Dur- 
ing the flare, this third component is hot (~1.5 keV), 
and cools to ~0.6 keV as the count rate decreases (see 
Figure Set 14 and Table 6). The flaring spectrum drives 
dramatic increases in both X-ray luminosity. Although 
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the flaring spectrum of Prox Cen exhibits a significant 
increase in the X-ray flux at wavelengths <10 A, the 
spectrum remains elevated (by an order of magnitude) 
compared to the quiescent spectrum out to long wave- 
lengths (see Figure 7). This is in contrast to earlier-type 
stars, which generally show a similar X-ray emission be- 
tween quiescent and flaring periods at wavelengths >100 


A. 


B. UPDATED X-RAY LUMINOSITIES FOR 
Swift/ XRT-IMAGED STARS 


We searched the Swift/XRT archive for all nearby 
stars that did not have available XMM-Newton or Chan- 
dra imaging. Nineteen of these stars (mostly K and M 
dwarfs; one system, Struve 2398, is a binary that is un- 
resolved in the Swift imaging) have fluxes or flux up- 
per limits available from Swift, which we summarize in 
Table 10. The R.A. and Decl. coordinates for each 
star is taken from Gaia Collaboration et al. (2023) and 
the stellar physical parameters are from Stassun et al. 
(2019). We use the quiescent spectrum from Table 6 
of the star that mostly closely matches each Swift star's 
stellar physical parameters and the Swift response files!! 
to convert Swift 0.3-10 keV count rates to luminosities. 
Figure 15 shows the Lx/Lyo ratio as a function of Teg 
(as in Figure 10) with the nineteen Swift stars added. 
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Table 10. Stars with Available Swift/XRT Imaging 


Star Alternate R.A. Dec. Spectral Distance Mass Teg = Lega logLx 
Name Name(s) (J2000) Type* (pc) (Mo) (K) (Le) ([ergs-i]) 
71 Ori HIP 29650 06:14:50.77 +19:09:20.29 F6V 20.87 1.85 6533 2.94 «21.59 
y Lep A GL 216A 05:44:27.45 -22:27:00.08 F6.5V 8.88 1.22 6258 2.34 27.53 
$? Cet GL 37 00:50:07.34 — -10:38:43.26 F7V 15.77 1.17 6250 1.85 <26.94 
7 CrB GL 584A 15:23:12.23 . +30:17:17.7 G2V 18.82 1.14 6029 1.64 «27.15 
61 Vir GL 506 13:18:24.97 -18:18:31.0 G7V 8.50 0.91] 5585 0.84 <26.78 
GL 567 LHS 5279  14:53:23.27 = +19:09:13.54 KOV 11.54 0.93 5258 0.54 28.28 
GL 68 LHS 1287  01:42:29.41 +20:15:55.87 K1V 7.60 0.88 5190 0.45 <26.29 
u Cas GL 53A 01:08:22.74 +54:54:47.53 K1V 7.55 0.82 5316 0.46 «26.71 
GL 785 LHS 488 . 20:15:18.88 = -27:02:01.61 K2V 8.79 0.85 5071 0.4 «26.50 
AK Lep GL 216B . 05:44:26.19 — -22:25:24.26 K2V 8.89 0.80 4869 0.3 28.00 
GL 688 HD 160346  17:39:16.72 +03:33:17.32 K2V 10.71 0.76 4982 0.34 «26.47 
Lacaille 8760 AX Mic 21:17:10.80 — -38:52:20.84 K9V 3.97 0.56 3599 0.1 26.89 
Wolf 1453 GL 205 05:31:28.21 — -03:41:11.50 MIV 5.70 0.56 3690 0.07 27.53 
Lacaille 9352 GL 887 23:06:00.94 — -35:50:49.79 MIV 3.29 0.54 3676 0.04 26.71 
BR Pis GL 908 23:49:13.59 +02:23:48.90 MIV 5.90 0.41 3685 0.029 <26.80 
GL 229A 0 LHS 1827  06:10:34.46  -21:52:04.16 — MI.5V 5.76 0.54 3912 0.06 27.55 
GL 1 HD 225213 00:05:31.99  -37:22:03.88 M1.5 4.36 0.38 3696 0.02 26.75 
Kruger 60 A GJ 860A 22:27:58.11  +57:41:38.52 M3V 4.01 0.34 3344 0.01 27.45 
Struve 2398B® GL 725B  18:42:43.94 +59:38:18.09 M3V 3.52 0.25 3345 0.021 27.57 
Struve 2398A% — GL 725A 18:42:43.94 +59:38:06.52 | M3.5V 3.52 0.33 3401 0.015 27.57 
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